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Abstract 


Multiple  input  multiple  output  (MIMO)  systems  can  provide  high  diversity, 
high  data  rate  or  a  mix  of  both,  for  wireless  communications.  This  report 
combines  both  modes  and  suggests  analyses  and  techniques  that  advance  the 
state  of  the  art  of  MIMO  systems.  Specifically,  this  report  studies  turbo 
space-time  coding  schemes  for  MIMO  systems. 

Before  the  designs  of  turbo  space-time  codes  are  presented,  a  fundamental 
tool  to  analyze  and  design  turbo  coding  schemes,  the  extrinsic  information 
transfer  (EXIT)  chart  method,  is  extended  from  the  binary/nonbinary  code 
case  to  coded  modulation  case.  This  extension  prepares  the  convergence 
analysis  for  turbo  space-time  code. 

Turbo  space-time  codes  with  symbols  precoded  by  randomly  chosen  uni¬ 
tary  time  variant  linear  transformations  (TVLT)  are  investigated  in  this  re¬ 
port.  It  is  shown  that  turbo  codes  with  TVLT  achieve  full  diversity  gain  and 
good  coding  gain  with  high  probability.  The  probability  that  these  design 
goals  are  not  met  is  shown  to  vanish  exponentially  with  the  Hamming  dis¬ 
tance  between  codewords  (number  of  different  columns).  Hence,  exhaustive 
tests  of  the  rank  and  the  determinant  criterion  are  not  required.  As  an  addi¬ 
tional  benefit  of  the  application  of  TVLT,  with  the  removal  of  the  constant 
modulation  condition,  it  is  proved  that  throughput  rates  achieved  by  these 
codes  are  significantly  higher  than  the  rates  achievable  by  conventional  space- 
time  codes.  Finally,  an  EXIT  chart  analysis  for  turbo  space-time  codes  with 
TVLT  is  developed,  with  application  to  predicting  frame  error  rate  (FER) 
performance  without  running  full  simulation. 

To  increase  the  data  rate  of  turbo-STC  without  exponentially  increas¬ 
ing  the  decoding  complexity,  a  multilevel  turbo  space-time  coding  scheme 
with  TVLT  is  proposed.  An  iterative  joint  demapping  and  decoding  receiver 
algorithm  is  also  proposed. 

For  MIMO  systems  with  a  large  number  of  transmit  antennas,  two  types 
of  layered  turbo  space-time  (LTST)  coding  schemes  are  studied.  For  systems 
with  low  order  modulation,  a  type  of  LTST  with  a  vertical  encoding  struc¬ 
ture  and  a  low  complexity  parallel  interference  cancellation  (PIC)  receiver 
is  shown  to  achieve  close  to  capacity  performance.  For  high  order  modu¬ 
lation,  another  type  of  LTST  with  a  horizontal  encoding  structure,  TVLT, 
and  an  ordered  successive  interference  cancellation  (OSIC)  receiver  is  shown 


to  achieve  better  performance  than  conventional  layered  space-time  coding 
schemes,  where  ordering  is  not  available  in  the  SIC  detection. 
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Chapter  1 
Introduction 


Wireless  communication  has  become  pervasive  in  every  day  life,  providing 
convenient  ways  to  exchange  information  as  well  as  supplying  emergency 
help  for  its  users  through  instant  accessibility.  Many  applications  such  as 
cellular  network,  wireless  LAN,  digital  video  broadcasting,  etc.,  have  been 
evolving  for  years  and  new  applications  are  still  emerging.  The  constant 
demand  for  wireless  communication  is  to  make  it  faster,  more  reliable,  and 
reaching  farther. 

One  of  the  biggest  obstacles  for  the  development  of  wireless  communica¬ 
tion  applications  is  the  fading  of  the  wireless  channel.  Fading  is  caused  by 
the  multipath  of  the  transmitted  signal  arriving  and  adding  either  in  phase  or 
out  of  phase  at  the  receiver,  leading  to  random  fluctuations  in  the  amplitude 
and  phase  of  the  received  signal. 

To  combat  fading,  multiple  diversity  technologies  such  as  frequency  di¬ 
versity,  temporal  diversity,  and  spatial  diversity,  have  been  proposed.  The 
core  idea  is  to  supply  the  receiver  with  multiple  copies  of  the  signal  so  that 
if  one  copy  of  the  signal  is  deeply  weakened  by  fading,  the  other  copies  may 
have  a  chance  to  be  strong. 

Among  the  diversity  techniques,  spatial  diversity  is  the  most  preferable 
since  it  does  not  cause  detection/decoding  delay  (as  time  diversity  does)  and 
does  not  cause  bandwidth  increase  (as  frequency  diversity  does).  Spatial 
diversity  is  achieved  by  employing  multiple  antennas  at  either  transmitter 
(transmit  diversity)  or  at  the  receiver  (receive  diversity),  or  both.  Wireless 
systems  with  multiple  transmit/receive  antennas  are  referred  to  as  multiple 
input  multiple  output  (MIMO)  systems. 

As  studies  advance,  it  is  shown  that  MIMO  systems  are  not  just  a  tech- 
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nique  to  overcome  the  problems  caused  by  fading.  In  fact,  the  studies  in 
information  theory  [2],  [3]  show  that  MIMO  systems  can  exploit  the  fading 
channels  and  achieve  much  higher  capacity  than  which  can  be  achieved  on 
a  nonfading,  additive  white  Gaussian  noise  (AWGN)  channel.  Space-time 
coding  (STC)  [1]  [4]  [5],  a  joint  design  of  coding,  modulation,  transmit  and 
receive  diversity,  has  been  shown  to  be  an  efficient  way  to  realize  the  big 
increase  in  the  spectrum  efficiency  predicted  by  the  information  theory. 

Turbo  codes  [6]  are  a  class  of  codes  that  achieves  near  capacity  per¬ 
formance.  They  have  also  been  extended  turbo  space-time  codes  (turbo- 
STC)  [7]  [8]  [9].  This  report  proposes  and  studies  a  new  type  of  turbo-STC 
with  the  time  varying  linear  transformation  (TVLT)  technique,  which  can 
be  designed  using  a  much  simpler  procedure,  and  achieves  a  higher  data  rate 
than  existing  turbo-STC. 

In  this  chapter,  the  capacity  of  a  MIMO  channel  is  first  presented,  fol¬ 
lowed  by  the  introduction  of  space-time  coding  and  turbo  coding.  Finally, 
turbo  space-time  coding  is  introduced  and  the  motivation  for  this  report  is 
presented. 


1.1  MIMO  Channel  Capacity 

Consider  a  MIMO  system  with  tit  transmit  and  Ur  receive  antennas.  The 
corresponding  channel  model  is 

r(f)  =  Hx(t)  +  n(t),  (1.1) 

where  r (t)  is  the  Ur  x  1  received  signal  vector  at  time  interval  t;  x(f)  is 
the  %  x  1  vector  of  the  transmitted  signal;  the  energy  of  of  each  symbols 
in  x(f)  is  Es/riT ;  H  =  {A,j}  is  the  ur  x  hr  channel  gain  matrix,  whose 
elements  hij , i  =  1,  j  =  1, ...,  n^,  are  complex  random  variables  due 

to  fading;  n (t)  is  zero-mean,  additive  white  Gaussian  noise  (AWGN)  with 
variance  Nq/2  per  dimension.  The  signal  to  noise  ratio  (SNR)  is  defined  as 
Es/N0. 

Consider  a  quasi  static  or  block  fading  channel  where  the  channel  gain 
matrix  H  keeps  constant  within  the  transmission  of  a  codeword  (frame)  but 
changes  independently  from  one  codeword  to  another  codeword.  The  capac¬ 
ity  of  the  MIMO  channel  (1.1)  for  a  realization  H  is  [2] 

C(H)  =  lcg2det(/  +  ;A-HHt) 
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Figure  1.1:  Example  of  linear  increasing  ergodic  capacity  with  number  of 
antennas. 

where  f  denotes  conjugate  transpose.  Since  H  is  random,  so  does  C  (H) .  The 
capacity  of  the  MIMO  fading  channel  can  be  described  as  ergodic  capacity 

C  (H)  p  (h1A, ...,  hnRtnT )  dh1A...dv. 

Consider  a  Rayleigh  fading  channel,  where  {hzj}  are  zero  mean  complex 
Gaussian  random  variables  with  variance  1/2  per  dimension.  Also  consider 
the  cases  where  nx  =  «r  =  n.  One  example  on  the  effect  of  the  number  of 
antennas  on  the  ergodic  channel  capacity  is  shown  in  Figure  1.1  for  SNR=10 
dB.  It  is  observed  that  the  channel  capacity  increases  linearly  with  the  num¬ 
ber  of  transmit  and  receive  antennas.  This  figure  demonstrates  the  advantage 
of  using  MIMO  systems  to  increase  the  data  rate. 

The  capacity  of  a  MIMO  fading  channel  can  also  be  described  by  the 
outage  probability,  defined  as  the  probability  that  the  channel  fails  to  support 
a  required  data  rate  Co 

Po  {Co)  =  Pr  (C  (H)  <  C0) 

Outage  probability  is  also  the  minimum  frame  error  rate  (FER)  that  can 
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Figure  1.2:  Outage  probability  comparison  of  channels  with  different  number 
of  transmit /receive  antennas. 


be  achieved  by  any  communication  system  over  a  certain  fading  channel. 
The  outage  probability  versus  SNR  for  Co  =  2  b/s/Hz  and  for  two  MIMO 
channels  with  nT  =  2,  Ur  —  1,  or  2  x  1,  and  nT  =  2,  ur  =  2,  or  2  x  2,  are 
shown  in  Figure  1.2,  with  a  comparison  to  the  the  outage  probability  of  the 
lxl  channel.  It  is  observed  that  the  slope  of  the  outage  probability  curve 
of  the  2x1  channel  is  twice  as  steep  as  that  of  the  lxl  channel,  showing 
the  effect  of  two  order  transmit  diversity.  The  slope  of  the  curve  for  the 
2x2  channel  is  steeper  than  that  of  the  2x1  channel,  showing  the  effect  of 
additional  receive  diversity. 


1.2  Space-time  Codes 

Encouraged  by  the  studies  in  information  theory,  space-time  coding  [1]  [4] 
[5]  has  been  developed  to  realize  the  data  rate  and  diversity  improvement 
predicted  by  those  studies.  Space- time  codes  are  multidimensional  codes  that 
generate  multiple  symbols  at  each  time  interval  to  be  transmitted  through 
multiple  antennas.  The  design  of  space- time  codes  is  guided  by  the  rank 
criterion  and  the  determinant  criterion.  Let  B  and  C  be  an  arbitrary  pair  of 
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Ut  x  L  codeword  matrices,  where  L  is  the  length  of  the  codeword.  Further 
let  E  =  B  —  C  be  the  error  matrix  between  B  and  C.  Then  for  the  Rayleigh 
fading  channel,  the  rank  criterion  says  that  the  space-time  code  achieves  r 
order  transmit  diversity,  if  and  only  if,  for  any  pair  of  codewords  B  and  C, 
the  rank  of  E  is  no  less  than  r.  The  determinant  criterion  says  that  if  r  order 
transmit  diversity  is  the  design  target,  to  achieve  maximum  coding  gain,  the 
minimum  value  of  the  product  of  the  determinants  of  the  all  r  x  r  principle 
cofactor  of  EE^  for  all  possible  pair  of  codewords,  should  be  maximized. 

These  two  criteria  reveal  the  fundamental  conditions  required  for  a  code 
to  achieve  a  desired  order  of  diversity  and  coding  gain.  However,  the  two 
criteria  are  based  on  a  pair  of  codewords.  That  is,  every  pair  of  codewords 
needs  to  be  tested  to  find  out  whether  the  two  criteria  are  satisfied.  This 
leads  to  difficulty  in  the  design  of  space-time  codes  and  a  computer  search  is 
usually  needed  [10]. 

The  rate  and  diversity  trade  off  is  also  studied  in  [1].  Let  tur  be  the 
diversity  advantage  of  the  system.  Assuming  that  the  signal  constellation  Q 
has  2b  elements,  the  rate  of  transmission  R  satisfies 

n  /  log  AbL(nT,r) 

R~ - L - 

in  bits  per  second  per  Hertz,  where  AbL  (nT,  r )  is  the  maximum  size  of  a  code 
length  L  and  minimum  Hamming  distance  r  defined  over  an  alphabet  of  size 
2bL.  In  a  special  case,  if  full  diversity  is  the  design  target,  then  the  maximum 
achievable  rate  is  b  b/s/Hz. 

As  an  example,  the  FER  performance  of  a  8  state  space-time  trellis  code 
(STTC)  [1]  with  4-PSK  modulation  and  rate  2  b/s/Hz  is  shown  in  Figure 
1.3,  with  comparison  to  the  outage  probability.  It  is  observed  that  both 
FER  curves  for  2  x  1  and  2x2  channels  are  parallel  to  those  of  the  outage 
probability,  demonstrating  that  full  diversity  is  achieved  by  this  STTC. 

The  decoding  complexity  for  STC  increases  exponentially  with  the  num¬ 
ber  of  transmit  antennas.  When  the  number  of  transmit  antenna  is  large 
(>  2),  layered  space-time  codes  (LSTC)  [5]  are  usually  used.  In  LSTC,  the 
transmit  antennas  are  divided  into  groups  with  smaller  number  (usually  two) 
antennas.  These  groups  are  also  referred  to  as  layers.  Independent  STC  are 
then  applied  in  each  layer.  At  the  receiver,  the  signals  from  different  layers 
are  first  separated  by  interference  suppression  techniques  and  then  decoded. 
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Figure  1.3:  FER  performance  of  a  8  state  STTC  [1]  with  comparison  to  the 
outage  probability. 


1.3  Turbo  Codes 

Turbo  code  [6]  is  a  strong  code  that  achieves  near  capacity  performance. 
It  is  constructed  by  either  parallel  or  serially  concatenating  two  component 
convolutional/block  codes.  An  example  of  a  parallel  concatenated  turbo 
code  is  shown  in  Figure  1.4.  In  the  encoder,  the  two  component  encoders 
are  connected  by  a  pseudo  random  interleaver  which  causes  the  generated 
codeword  to  be  pseudo  random  also.  This  is  the  reason  why  turbo  codes  can 
achieve  near  capacity  performance.  Usually,  the  interleaver  should  be  long 
enough  for  the  turbo  code  to  achieve  the  superior  performance. 

After  being  encoded  by  the  two  component  encoders,  the  information 
bits  and  the  parity  bits  may  be  mapped  into  nonbinary  symbols.  These 
types  of  schemes  are  referred  to  as  turbo  coded  modulation  [11]  [12].  If  the 
multiple  symbols  are  generated  after  the  mapping  and  transmitted  through 
multiple  antennas,  the  turbo  code  will  be  referred  to  as  turbo  space-time 
codes  [7]  [8]  [9]. 

The  decoder  of  the  turbo  code  is  also  composed  of  two  component  soft 
input  soft  output  (SISO)  decoders,  and  works  in  an  iterative  fashion.  Dur¬ 
ing  the  iterations,  the  performance  of  the  two  component  (SISO)  decoders  is 
enhanced  by  exchanging  with  each  other  the  extrinsic  information,  or  newly 
generated  information,  about  the  transmitted  data.  The  convergence  per¬ 
formance  of  the  iterative  decoder  can  be  conveniently  analyzed  using  the 
extrinsic  information  transfer  chart  (EXIT)  [13]  and  other  similar  meth¬ 
ods  [14],  [15]. 

1.4  Turbo  Space-time  Codes 

Turbo  space-time  codes  (turbo-STC)  [7]  [8]  [9]  are  a  special  type  of  turbo 
codes  that  are  applied  to  MIMO  systems.  As  described  in  the  last  section, 
each  pair  of  the  codewords  should  be  tested  to  see  if  the  rank  criteria  and 
determinant  criterion  are  satisfied.  The  test  is,  however,  even  more  difficult 
to  implement  than  the  test  for  STTC.  Since  usually  turbo  space-time  codes 
have  a  long  interleaver  (codeword),  the  number  of  codewords  to  test  will  be 
very,  if  not  prohibitively,  large.  This  is  one  of  the  problems  to  be  dealt  with 
in  this  report.  A  technique  called  time  varying  linear  transformation  (TVLT) 
will  be  applied  to  guarantee  the  full  diversity  of  turbo-STC  without  testing 
each  pair  of  the  codewords. 


12 


The  EXIT  chart  method  is  an  accurate  and  low  complexity  method  used 
to  predict  the  performance  of  turbo  codes.  But  it  is  not  useful  for  conven¬ 
tional  turbo-STC  for  the  reasons  that  will  be  explained  later  in  section  3.3.4 
of  this  report.  With  TVLT,  this  method  can  be  again  applied.  Generating 
EXIT  charts  for  turbo  space-time  codes  with  TVLT  is  another  problem  to 
be  investigated  in  this  report. 


1.5  Outline  of  the  report 

This  report  is  organized  as  follows.  First  the  EXIT  chart  is  extended  from 
binary  turbo  codes  to  turbo  coded  modulation  in  Chapter  2,  which  prepares 
the  tool  for  analyzing  the  convergence  of  turbo  space-time  codes.  In  Chapter 
3,  the  technique  of  TVLT  is  applied  to  turbo-STC.  It  will  be  shown  that 
TVLT  guarantees  full  diversity,  and  the  large  minimum  Hamming  distance 
of  turbo  space-time  codes  will  guarantee  good  coding  gain.  The  convergence 
analysis  of  turbo-STC  is  also  presented  in  Chapter  3.  Multilevel  turbo-STC 
and  layered  turbo-STC  will  be  presented  in  Chapters  4  and  5,  respectively. 
The  summary  and  anticipated  future  work  are  presented  in  Chapter  6. 
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Chapter  2 

Convergence  Analysis  of  Turbo 
Trellis  Coded  Modulation 


Convergence  analysis  of  iterative  decoding  algorithms  for  turbo  codes  has 
received  much  attention  recently  due  to  its  useful  application  in  predicting 
code  performance,  its  ability  to  provide  insights  into  the  encoder  structure, 
and  its  usefulness  in  helping  with  the  code  design.  Several  methods  have 
been  proposed  to  analyze  the  convergence  of  the  iterative  decoders  for  bi¬ 
nary  turbo  codes  [13-15].  In  particular,  the  extrinsic  information  transfer 
(EXIT)  method  [13]  has  created  a  lot  of  interest.  The  EXIT  method  has 
been  extended  to  the  analysis  of  nonbinary  turbo  codes  in  [16]. 

Turbo  trellis  coded  modulation  (TTCM)  [11,12],  conjoins  signal  map¬ 
ping  techniques,  such  as  Ungerboeck’s  signal  space  partition  [17],  with  turbo 
coding,  to  achieve  significant  coding  gains  without  increasing  bandwidth. 
However  the  need  for  signal  mapping  makes  the  encoder  structure  more  com¬ 
plex  to  design  and  analyze  than  binary  turbo  codes.  Hence  the  convergence 
analysis  is  a  very  important  tool  for  the  design  and  comparison  of  TTCM 
schemes. 

In  TTCM,  systematic  bits  and  parity  bits  are  usually  mapped  to  and 
transmitted  as  a  single  symbol.  At  the  receiver,  a  symbol  decoder  is  used 
instead  of  a  binary  decoder  [11, 12] .  Consequently,  the  systematic  information 
and  extrinsic  information  are  not  naturally  separated  as  in  binary  turbo 
codes.  The  separation  is  not  necessary  for  carrying  out  the  decoding.  It 
is,  however,  needed  for  generating  the  EXIT  chart.  Due  to  mixing  of  the 
systematic  and  extrinsic  information,  the  approach  for  generating  the  EXIT 
chart  for  binary  decoders,  or  even  for  nonbinary  decoders  [16],  can  not  be 
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directly  applied  to  TTCM. 

In  this  chapter,  the  EXIT  chart  for  TTCM  is  developed  by  showing  a  way 
to  explicitly  separate  the  systematic  information  from  the  extrinsic  informa¬ 
tion.  To  be  specific,  the  analysis  based  on  the  TTCM  codes  introduced  in  [11] 
is  developed.  The  application  of  the  proposed  method  to  other  structures 
is  fairly  straightforward.  Furthermore,  the  EXIT  chart  method  presented 
in  this  chapter  will  be  used  to  analyze  the  convergence  performance  of  the 
turbo  space-time  codes. 

The  rest  of  the  chapter  is  organized  as  follows.  In  section  2.1  the  structure 
of  TTCM  is  briefly  introduced  and  the  EXIT  chart  for  TTCM  is  developed 
in  section  2.2.  Several  examples  are  used  to  illustrate  and  demonstrate  the 
proposed  method  in  section  2.3.  Conclusions  are  drawn  in  section  2.4. 


2.1  TTCM  Structure 

2.1.1  Encoder  Structure 

The  structure  of  the  TTCM  encoder  [11]  is  presented  in  Figure  2.1.  Assuming 
a  2m+1  constellation,  at  each  time  interval  k,  a  symbol  is  transmitted  repre¬ 
senting  m  information  bits  b [k]  =  {b[k,  i]}™^  •  Of  the  m  bits,  possibly  only 
fh  bits  {b[k,  f]}™g\  are  encoded,  whereas  the  other  (m  —  fh)  uncoded  bits  are 
input  directly  to  the  symbol  mapper.  The  turbo  encoder  consists  of  two  rate 
fh/(fh+  1)  component  recursive  systematic  convolutional  (RSC)  encoders, 

which  are  connected  by  a  symbol  (fh  bits)  interleaver.  These  fh  encoded 

fh 

bits,  are  labeled  for  later  use  as  the  2m-ary  level  symbol  d[k\  =  EMM]* 

z=l 

The  parity  bits  C\[k)  and  C2 [A]  from  the  two  component  encoders  (the  parity 
bits  from  the  encoder  with  interleaved  input  are  deinterleaved  to  their  orig¬ 
inal  order)  are  alternatively  punctured,  and  the  surviving  parity  bit  c[k]  is 
mapped  together  with  the  m  data  bits  to  M-QAM  or  MPSK  symbols  s[k], 
and  transmitted. 
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b[k,m-\] 

b[k,m\ 

b[k,m-\] 
b[k,  0] 


Figure  2.1:  Structure  of  TTCM  encoder  with  m  coded  bits  and  m  —  rh 
uncoded  bits. 

2.1.2  Decoder  Structure 

The  structure  of  the  TTCM  decoder  is  depicted  in  Figure  2.2.  The  decoder 
consists  of  two  symbol  maximum  a  posteriori  (MAP)  decoders,  which  ex¬ 
change  information  about  the  symbols  {d[k]}  formed  by  the  encoded  data 
bits.  Of  course,  no  information  on  uncoded  bits  is  exchanged.  Hence  in 
TTCM,  the  systematic  information  represents  the  rh  coded  bits  {b[k,i} ■ 
Let  the  TTCM  codeword  consist  of  N  symbols  (s[l],  s[2], . . . ,  s[N]).  The 
channel  observation  is  given  by 

y[k]  =  s[k]  +  n[k],  (2.1) 

where  n[k ]  is  zero  mean  white  Gaussian  noise  with  variance  a2.  The  signal 
to  noise  ratio  (SNR)  is  then  defined  as  E[\s[k]\2]/ a2 . 

Due  to  fact  that  the  symbol  s[k]  embodies  both  the  systematic  and  par¬ 
ity  bits,  and  consequent  to  the  alternative  puncturing  of  the  parity  bits,  it  is 
shown  in  [11]  that  half  of  the  information  exchanged  between  component  de¬ 
coders  contains  systematic  information  and  half  does  not.  Unlike  the  binary 
turbo  code  case,  the  a  posteriori  probabilities  (APP)  Lp[fc]  =  {log  P(d[k]  — 

0m 

■u|y)}u=Q  computed  by  either  component  decoder  can  be  decomposed  into 
only  two  parts.  The  first  part  is  the  a  priori  input  La[A:]  A  {log  P(d[k]  = 

u)}w=o  1>  passed  from  the  other  decoder.  The  second  part  is  different  for  two 
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cases:  (i)  at  time  interval  k ,  s[k]  contains  the  parity  bit  from  the  correspond¬ 
ing  component  encoder;  (ii)  at  time  interval  k ,  s[k]  contains  parity  bit  from 
the  other  component  encoder.  For  case  (ii),  the  second  part  only  contains 
the  newly  generated  extrinsic  information  output,  denoted  as  L°ut  [/c] .  For 
case  (i)  however,  the  second  part  is  a  mixture  of  the  extrinsic  information 
and  systematic  information,  denoted  as  L “*[&].  That  is, 


/  La[fc] +L2?[fc],  k  E  case  (i) 

\  L0[fc]  +  L °ut[k],  k  e  case  (ii)  ’ 


(2.2) 


Accordingly,  the  a  priori  input  to  each  component  decoder  La[/c]  is,  for  case 
(i),  only  extrinsic  information,  denoted  as  L” n[k]\  and  is,  for  case  (ii),  the 
mixture  of  both  extrinsic  information  and  systematic  information,  denoted 
as  L’"[/c].  That  is 


f  L*n  [k] ,  k  e  case  (i) 
1  k  e  case  (“) 


(2.3) 


L.  =  L  or  L 

X 


Dec  1 


i  Ii  ~  I.  +L.,  or  I.  + 1. 
Symbol  MAP  "  “  "  ‘ 


L„  =  L„  +  L_  or  La  +  L 

p  a  es  a  € 


1  Symbol  MAP 

A 

h= 

-0  + 

I  IL  | 


a  e  «s 


Figure  2.2:  The  structure  of  a  TTCM  decoder  with  two  component  MAP 
decoders. 


2.2  EXIT  Charts  for  TTCM 

The  EXIT  chart  is  a  tool  for  studying  the  convergence  of  turbo  decoders 
without  simulating  the  whole  decoding  process.  The  chart  is  generated  by  a 
single  component  decoder  by  mimicing  the  a  priori  input  (as  if  it  is  generated 
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by  the  other  decoder)  and  measuring  the  output.  It  is  found  in  [13-15]  that 
for  binary  MAP  decoders,  when  the  random  interleaving  sufficiently  deep,  the 
a  priori  inputs  (the  extrinsic  information  passed  from  the  other  component 
decoder),  A  =  log  ,  for  each  data  bit  6,  are  mutually  independent  and 
can  be  closely  modeled  as  independent,  identically  distributed  (i.i.d.)  random 
variables  that  follow  the  same  Gaussian  distribution  p(A|fe).  In  addition,  the 
mean  value  p  and  variance  er2  of  p(A|6)  satisfy  the  consistency  condition 
p  =  \cr2.  Then  p(A|6)  is  determined  by  only  one  parameter,  either  p  or  a2. 
So  to  generate  A  with  a  given  mutual  information  content  A  =  the 

parameter  of  p(A|6),  p  (or  a2)  can  be  found  by 


A  = 


A 


/(A;  b) 

1  _  r+oo 

-J2  P(\\b=x) 

Z  i=0,l  •'_0° 

,log  2-p(A|6  =  x) 

W 


p(X\b  =  0)  +  p(A|5  =  1) 


d\ 


(2.4) 


and 

p  =  J-\A ).  (2.5) 

The  above  procedure  is  used  to  generate  the  a  priori  input  to  a  binary 
MAP  decoder.  For  TTCM,  the  component  decoder  is  a  nonbinary  symbol 
decoder,  and  the  a  priori  has  different  contents  in  actual  operation  for  the 
two  different  cases  shown  in  (2.3).  To  mimic  the  a  priori  input  UJ1 [k]  for 
k  €  case  (i),  which  contains  only  the  extrinsic  information  generated  by 
the  other  decoder,  the  assumption  proposed  in  [16]  that  the  log-likelihood 
ratios  (A[A,  i]  =  log  P(fc[fc|Sb)  }£=oX  f°r  the  bits  {b[k,  zj}’"o1  composing  the 
nonbinary  symbol  d[k]  are  i.i.d.  Gaussian  random  variables,  is  adopted. 
Thus  { A[A,  i]  can  he  generated  as  in  the  binary  turbo  code  case,  with 

the  value  A  in  (2.5)  set  to  A  =  I^/rh,  where  I™  is  the  desired  mutual 
information  between  L ™\k\  and  d[k],  i.e.,  /*"  =  I(L™[k)\d[k]).  The  term 

L ’"[&]  —  (log  P(d[k]  =  m)}^^1  can  be  calculated  using  {AfA:,^}^1. 

For  k  €  case  (ii),  the  a  priori  is  a  mixture  of  both  extrinsic  and  system¬ 
atic  information  L*"[fc]  as  shown  in  (2.3).  Before  generating  the  mixture,  it 
is  needed  to  find  out  what  the  systematic  information  will  be,  even  if  does 
not  appear  explicitly  during  the  decoding  process  in  TTCM.  In  binary  turbo 
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coding,  the  systematic  information  is  the  information  associated  with  the 
channel  observation  of  the  systematic  bit.  This  should  also  be  the  case  in 
TTCM,  although  each  sample  of  the  channel  observation  contains  both  sys¬ 
tematic  bits  (symbol),  uncoded  bits  and  parity  bits.  That  is,  the  systematic 
information  about  symbol  d[k\  should  be  computed  as 

pm  =  *i#])=  E  -  E  E  '’«#[*!  = 

6[fc,m+l]={0,l}  bk>m={ 0,1}  cfc={0,l} 

mapping(d[/c]  =  u,  b[k,  m  +  1],  ...,&[&,  ra],c[A;])),  u  =  0, ...,  2m  —  1.  (2.6) 

Equation  (2.6)  can  also  be  rewritten  in  the  form  of  a  vector  of  log- 
likelihood  values  as 


U[k]  =  {log P(d[k]  =  u\y[k])yu  =  0, ...,  2*  -  1}.  (2.7) 

Now  the  systematic  information  Ls[fc],  for  k  G  case  (ii),  has  been  explicitly 
formed.  To  generate  the  mixture,  it  is  also  necessary  to  generate  for  each 
k  G  case  (ii)  a  vector  extrinsic  part  L *n[fc]  with  /(L*n [/,:];  d\k])  —  7*n,  in  the 
same  way  as  for  case  (i).  Thus,  the  mixture  L0[fc]  =  L *"[&]  is  generated  as  the 
combination  of  extrinsic  information  and  systematic  information,  i.e.,  L a[k\  ~ 
L *n[/c]  +  Ls[/c].  To  conclude,  the  a  priori  input  to  the  TTCM  component 
decoder  is  generated  as 

L*n[/c],  k  6  case  (i) 

L ‘"[A:]  +Ls[fc],  k  G  case  (ii). 

Next  the  measurement  of  the  extrinsic  information  at  the  output  of  the 
TTCM  component  decoder  is  discussed.  According  to  (2.2),  for  k  G  case 
(i)  the  extrinsic  information  is  mixed  with  systematic  information.  However, 
the  systematic  information  can  be  formed  using  (2.6)  and  (2.7)  and  then 
subtracted.  That  is,  the  extrinsic  information  L°ut  [k]  for  data  symbol  d[k]  is 
obtained  as: 


f  L p[k]  -  La[fc]  -  Ls[k],  k  G  case(i) 

\  Lp[/c]  —  La[fc],  k  G  case  (ii). 


(2.9) 


Due  to  the  independence  assumption,  {L^fA;]}^1  are  mutually  inde¬ 
pendent  and  follow  the  same  distribution  pe  (L^  |  d) .  The  multi-dimensional 
joint  distribution  pe(L°u(|d)  is  obtained  by  measuring  the  histogram  of 
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using  its  samples  {L^fc]}^1.  Then  the  mutual  information  between  the 
extrinsic  output  and  the  transmitted  symbol  d  is  computed  as 


rout 


log 


(2ft-f)-Prrf(g[rf  =  u)JC 
E^o  lpTim=u) 


(2.10) 


The  EXIT  chart  for  TTCM  is  then  generated  by  measuring  the  output 
mutual  information  I°ut  as  a  function  of  a  family  of  inputs  with  specified 
mutual  information  Pen  and  channel  observation  of  specified  SNR,  I°ut  = 
/(/f,  SNR). 

The  overall  block  diagram  for  generating  EXIT  charts  for  TTCM  schemes 
is  illustrated  in  Figure  2.3. 
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Figure  2.3:  The  block  diagram  to  generate  EXIT  charts  for  a  TTCM  MAP 
decoder. 


2.3  Applications 

In  this  section,  a  few  examples  are  presented  to  illustrate  the  applications  of 
the  proposed  EXIT  chart  method  for  TTCM. 
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2.3.1  Role  of  Primitive  Feedback  Polynomials 

The  proposed  EXIT  chart  method  is  applied  to  analyze  and  compare  two 
TTCM  schemes.  The  first  code  is  from  [11],  Table  I,  first  row.  This  is  an 
8-state  code  with  8-PSK  modulation,  and  m  =  in  =  2.  The  second  code  has 
the  same  structure  but  the  parameter  H°(D )  in  the  first  row  of  table  I  [11] 
is  changed  from  (11)8  to  (13)g.  This  changes  the  feedback  polynomial  of  the 
component  encoder  from  nonprimitive  to  primitive.  Hence,  a  primitive  code 
and  a  nonprimitive  code  are  compared.  The  data  rates  of  both  codes  are  2 
b/s/Hz. 

The  EXIT  charts  for  these  two  codes  are  presented  in  Figure  2.4,  for 
channel  SNR  5.8  dB,  6.2  dB  and  7.0  dB,  respectively.  For  schemes  with  sym¬ 
metric  structure  (identical  component  encoders/decoders),  only  the  EXIT 
chart  for  one  of  the  decoders  needs  to  be  plotted.  In  the  figure,  the  abscissa 
is  PJ1  =  /(L™;  d)  in  bits,  and  the  ordinate  is  I°ut  —  I (L°ut ;  d )  also  in  bits. 

It  is  observed  from  the  figure  that  for  both  codes,  the  curves  of  SNR=5.8 
dB  intersect  with  the  y  =  x  line,  which  means  that  at  some  point  during 
the  iteration  the  decoder  can  no  longer  increase  the  extrinsic  information 
compared  to  the  input.  In  this  case  the  decoding  will  fail  no  matter  how 
many  iterations  are  run.  Both  curves  for  SNR  =  6.2  dB  are  just  a  little 
above  the  y  =  x  line,  which  means  for  that  any  SNR  ^  6.2  dB,  the  extrinsic 
information  output  at  the  component  decoder  will  always  be  better  than 
the  input,  and  the  decoder  will  converge  to  the  correct  codeword  with  high 
probability.  That  is,  the  convergence  threshold  of  these  two  codes  is  about 
6.2  dB.  The  simulated  bit  error  rate  (BER)  performance  curve  is  shown  in 
Figure  2.5  after  8  iterations.  The  length  of  the  interleaver  is  1024.  It  is 
observed  that  the  BER  curve  actually  turns  down  sharply  at  about  6.2  dB, 
which  verifies  the  convergence  threshold  predicted  by  the  EXIT  chart. 

In  Figure  2.4,  it  is  also  observed  that  the  curves  of  the  primitive  code 
are  always  higher  than  the  nonprimitive  code,  especially  at  higher  values  of 
input  mutual  information.  This  implies  a  faster  convergence  speed,  which 
agrees  with  the  observation  in  binary  turbo  codes  that  codes  with  primitive 
feedback  polynomial  have  faster  convergence  speed  than  the  nonprimitive 
ones  [15].  This  observation  is  also  verified  by  the  BER  curves  in  Figure  2.5, 
where  with  the  same  number  of  iterations,  the  performance  of  the  primitive 
code  is  always  better  than  that  of  the  nonprimitive  code.  So  this  example 
suggests  that  in  TTCM,  primitive  codes  should  be  favored  over  nonprimitive 
codes. 

I 
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Figure  2.5:  BER  comparison  of  two  codes  with  primitive  and  nonprimitive 
feedback  polynomials. 

2.3.2  Analysis  and  design  of  TTCM  schemes  with  un¬ 
coded  bits 

In  this  subsection,  the  method  to  use  the  proposed  EXIT  scheme  to  predict 
the  performance  of  TTCM  schemes  with  uncoded  bits  is  shown.  It  is  also 
shown  how  the  system  design  can  be  improved  by  examining  the  EXIT  chart. 

First  the  3  b/s/Hz,  16  QAM  scheme  in  [11],  with  8-state  component  codes 
is  studied.  In  this  scheme  three  bits  are  transmitted  at  each  time  interval 
but  only  two  are  coded.  Consequently,  there  are  two  parallel  branches  for 
each  transition  in  the  trellis.  The  goal  of  the  analysis  is  to  find  at  what  SNR 
the  system  can  achieve  a  given  bit  error  rate  (BER),  e.g.,  BER=10~5.  First 
according  to  the  set  partition  diagram,  [17]  Figure  5,  the  distance  between 
the  parallel  transitions  is  found  to  be  dp  =  \/8  ■  A0,  where  A0  is  the  minimum 
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distance  between  the  points  of  16  QAM  and  A0  =  2/VTO  •  Es.  Since 


BER 


the  required  SNR  is 


N0 


=  ^  (Q~'  (BER))2 


(2.11) 


(2.12) 


For  BER=10-5,  the  required  channel  SNR  is  about  10.5  dB.  Next,  the  EXIT 
chart  for  the  underlying  turbo  code  at  this  SNR  is  generated  and  shown 
in  Figure  2.7.  From  the  figure  it  is  observed  that  EXIT  curve  is  far  above 
the  y=x  line  and  the  opening  between  the  EXIT  curve  and  the  y=x  line  is 
very  large,  which  means  the  decoder  has  converged  at  SNR=10.5  dB  and  the 
convergence  speed  is  fast.  The  convergence  threshold  for  this  code  is  found 
as  9.6  dB,  as  also  shown  in  Figure  2.7.  Obviously,  the  BER  performance 
of  the  whole  system  is  dominated  by  the  BER  of  the  uncoded  bit  and  the 
performance  of  BER=10-5  will  be  achieved  at  about  10.6  dB.  This  can  be 
verified  by  checking  the  BER  performance  curve  in  [11]. 

The  analysis  also  suggests  the  turbo  code  used  is  too  strong,  and  a  code 
with  lower  state  complexity  may  be  used.  In  Figure  2.6  A  4-state  component 
code  is  proposed  as  an  alternative.  The  EXIT  chart  for  this  code  is  also 
shown  in  Figure  2.7,  for  SNR=10.5  dB.  It  is  observed  that  although  the 
curve  for  the  4-state  code  is  lower  than  the  8-state  code,  the  opening  is  still 
very  wide  and  suggests  a  good  convergence  speed.  Hence,  the  complexity  of 
the  decoder  can  be  lowered  with  little  performance  loss. 


2.3.3  Performance  Analysis  of  TTCM  over  Fading  Chan¬ 
nel 

The  convergence  analysis  of  the  TTCM  schemes  over  AWGN  channel  also 
leads  to  convenient,  accurate  prediction  of  the  performance  of  the  TTCM 
schemes  over  block  fading  channels.  The  frame  error  rate  (FER)  perfor¬ 
mance  is  just  the  probability  that  the  channel  gain  falls  below  the  conver¬ 
gence  threshold  of  the  TTCM  code.  For  the  block  Rayleigh  fading  channel, 
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Figure  2.6:  An  alternative  4-state  component  encoder  for  the  3  b/s/Hz, 
16QAM  scheme. 


Figure  2.7:  EXIT  charts  for  two  16  QAM,  3  b/s/Hz  TTCM  schemes  with 
8-state  and  4-state  component  encoder,  respectively. 


the  channel  model  is 


y[k\  =  h  ■  s[k] +n[k\,  (2.13) 

where  h  is  the  complex  channel  gain  modelled  as  a  zero  mean,  unit  variance 
Gaussian  random  variable.  The  definition  of  other  variables  are  the  same 
as  in  (2.1).  Denote  the  channel  SNR  7  =  The  effective  channel  SNR 
under  channel  gain  h  is  £  =  7  ■  \h\2.  It  is  obvious  f  follows  an  exponential 
distribution.  The  probability  density  function  of  7  is 

p(0  =  ^  exp(-£/7)  (2-14) 

Denote  the  SNR  required  by  the  TTCM  code  to  converge  -yth.  Then  the  FER 
is 


Pe  =  P(  7<7tfc)  (2.15) 

Pth  1 

=  ~  e*p(— £/7K 

Jo  7 

=  1  —  exp(— 7^/7) 

The  convergence  threshold  for  the  primitive  2  b/s/Hz  TCM  code  dis¬ 
cussed  in  subsection  2.3.1  is  7t/l  =  6.2  dB.  The  predicted  FER  performance 
versus  SNR  is  plotted  in  Figure  2.8.  The  FER  performance  obtained  by  full 
simulation  is  also  shown  in  the  figure  for  comparison.  It  is  observed  that  the 
two  curves  almost  overlap  each  other,  indicating  an  accurate  match. 


2.4  Chapter  Summary 

The  method  to  generate  the  EXIT  chart  for  TTCM  decoders  is  proposed. 
The  accuracy  and  usefulness  of  this  method  has  been  verified  through  exam¬ 
ples.  Since  the  complexity  of  generating  EXIT  charts  is  significantly  lower 
than  a  full  simulation,  the  proposed  method  supplies  a  convenient  way  to 
compare  various  TTCM  schemes  as  well  as  to  verify  and  compare  new  de¬ 
signs. 
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Comparison  of  the  simulated  and  predicted  FER  performance  for 
TCM  code. 
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Chapter  3 

Turbo  Space-time  Coding  with 
Time  Varying  Linear 
Transformat  ions 


The  invention  of  space-time  codes  [1]  has  been  seminal  in  the  advance  of 
multiple  input  multiple  output  (MIMO)  techniques.  Space-time  codes  apply 
to  multiple  transmit  antennas  to  achieve  both  diversity  and  coding  gains 
over  fading  channels.  However,  the  design  of  full  diversity  codes  with  good 
coding  gain  has  been  challenging.  One  of  the  difficulties  has  been  the  lack  of 
effective  design  tools  to  meet  specified  diversity  and  coding  gains  criteria.  In 
most  cases,  to  ensure  meeting  diversity  and  gain  specifications  [1],  exhaustive 
checks  are  required  of  all  pairs  of  codewords,  and  numerical  simulations  are 
needed  to  evaluate  overall  performance. 

Turbo  space-time  codes  (STC)  are  specially  designed  turbo  codes  for 
MIMO  channels  that  can  achieve  much  higher  coding  gains  than  conventional 
space-time  trellis  codes.  Existing  schemes  such  as  [7],  [8],  and  [9],  etc.,  display 
similar  performance  but  have  distinct  encoder/decoder  structures.  For  all 
turbo-STC,  exhaustive  tests  are  required  to  ensure  full  diversity.  If  the  test 
fails,  the  codes  have  to  be  adjusted,  and  the  testing  procedure  is  restarted. 

The  main  goal  in  this  chapter  is  to  simplify  the  design  of  full  diversity 
turbo  space-time  codes.  Phase  sweeping  [18-20],  a  precursor  to  space-time 
codes,  guarantees  full  diversity,  but  no  coding  gain.  It  still  attracts  attention 
both  in  academia  and  industry  [21,22] .  The  approach  taken  in  this  report  is  a 
technique  related  to  phase  sweeping  and  is  referred  to  as  time  varying  linear 
transformation  (TVLT)  [23].  TVLT  offers  more  flexibility  for  integration  in 
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space-time  codes  than  phase  sweeping.  With  TVLT,  transmitted  symbols  are 
modified  by  a  unitary  transformation  represented  by  a  unitary  matrix.  The 
elements  of  the  matrix  are  in  part  selected  at  random  and  in  part  designed 
to  meet  the  unitary  conditions.  The  transformation  changes  at  each  time 
interval,  but  the  sequence  of  transformations  is  known  to  both  the  transmitter 
and  the  receiver.  A  previous  attempt  of  applying  TVLT  to  turbo  space-time 
codes  proved  unsuccessful. 

It  will  be  shown  that  the  probability  of  losing  full  diversity  and  not  meet¬ 
ing  the  coding  gain  specification  vanishes  exponentially  with  the  Hamming 
distance  between  codewords  defined  as  the  number  of  different  columns  in 
pairs  of  codewords.  Since  turbo  codes  typically  feature  considerable  Ham¬ 
ming  distances,  an  exhaustive  check  of  all  pairs  of  codewords  is  not  necessary. 

Besides  simplifying  design,  TVLT  can  also  be  used  to  construct  full  diver¬ 
sity  space-time  codes  with  higher  rate  than  codes  with  constant  modulation, 
e.g.,  4-PSK,  at  each  time  interval.  The  maximum  data  rate  limit  derived 
in  [1]  for  a  full  diversity  space-time  codes  is  based  on  the  assumption  of  con¬ 
stant  modulation.  TVLT  has  the  effect  of  varying  the  modulation  at  each 
time  interval.  It  leads  to  a  higher  maximum  rate  that  will  be  derived  in  this 
report.  For  example,  full  diversity  space-time  codes  with  8-PSK  modulation 
have  a  maximum  rate  of  3  b/s/Hz.  A  full  diversity  TVLT  8-PSK  code  with 
a  rate  of  4  b/s/Hz  will  be  demonstrated. 

Finally,  it  is  shown  that  turbo  STC  with  TVLT  lend  themselves  to  the 
application  of  a  form  of  the  extrinsic  information  transfer  (EXIT)  chart 
especially  modified  for  the  task.  This  method  is  an  extension  of  the  EXIT 
chart  for  single  antenna  turbo  codes  [13, 16, 24]  to  turbo  STC.  It  provides 
an  efficient,  low  complexity  method  to  quickly  predict  the  performance  of 
turbo  STC.  With  this  convenient  tool,  different  designs  and  different  codes 
and  can  be  easily  compared,  and  the  best  one  can  be  selected.  Several  new 
turbo  STC  are  developed  and  analyzed  in  this  manner. 

The  remainder  of  this  chapter  is  organized  as  follows.  In  Section  3.1,  after 
presenting  the  channel  model,  TVLT  codes  are  introduced.  Salient  properties 
of  turbo  STC  with  TVLT  are  explored  in  Section  3.2.  The  performance 
analysis  based  on  the  EXIT  charts  is  developed  in  Section  3.3.  Numerical 
results  are  presented  in  Section  3.4  and  concluding  remarks  are  found  in 
Section  3.5. 
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3.1 


Channel  Model  and  TVLT 


In  this  chapter  MIMO  systems  with  nT  transmit  and  nR  receive  antennas 
over  block  Rayleigh  fading  channel  are  considered.  The  received  signal  at 
time  interval  t  is 

r  (t)  =  Hx(f)  +  n  (t),  (3.1) 

where  x(t)  is  the  nT  x  1  vector  of  the  transmitted  signal;  r (t)  is  the  nR  x  1 
received  signal  vector;  H  is  the  nR  x  nr  channel  gain  matrix,  which  con¬ 
sists  of  complex- valued  scalars  hij,i  =  1  ,...,nR,  j  =  1  ,...,nT,  modeled  as 
zero-mean,  mutually  independent,  identically  distributed  Gaussian  random 
variables  with  unity  variance  The  term  n  (t)  is  zero-mean,  additive  white 
Gaussian  noise  (AWGN)  with  variance  nT/  (2  SNR)  per  dimension,  where 
SNR  is  the  signal  to  noise  ratio  per  symbol.  With  this  model,  the  mean  SNR 
at  the  receiver  is  independent  of  the  number  of  transmit  antennas. 

The  transmitter  with  TVLT  is  shown  in  Figure  3.1.  Denote  a  codeword 
produced  by  the  encoder  C  =  [c(l), ...,  c(L)],  with  L  the  codeword  length. 
At  each  time  interval  t,t  =  1, ..,  L,  a  vector  of  coded  symbols  c(t)  —  [ci(f), ..., 
CnT(t)}T  is  multiplied  by  a  predetermined  nT  x  nr  unitary  matrix  G(t)  and 
yields  the  vector  of  symbols  x(t),  which  is  subsequently  transmitted, 

x(t)  =  G(t)c(t).  (3.2) 

The  set  of  unitary  matrices  {G(t)}^=1  is  generated  randomly,  independent 
of  the  codewords  it  is  applied  to.  Once  generated,  the  set  is  stored  both 
at  the  transmitter  and  at  the  receiver  and  it  stays  fixed  for  the  duration  of 
the  communication.  Insofar  as  the  elements  of  G(t)  are  time- varying  and 
realizations  of  random  variables,  x(t)  appear  to  be  subject  to  a  time-varying, 
random  modulation.  At  the  receiver,  a  random  modulation  is  considerably 
simpler  to  demodulate  than  to  decode  a  random  code  since  demapping  is 
carried  out  symbol  by  symbol  rather  than  on  long  sequences. 

Next,  the  generation  of  these  unitary  matrices  for  an  arbitrary  number  of 
transmit  antennas  is  discussed.  For  p,  q  =  1, ...,  nr,  let  gp,q{t)  be  the  elements, 
and  g g(t)  be  the  columns  of  G(f).  Since  G(t)  is  unitary,  it  must  satisfy 

4(*)*.w-{i  .  <3.3) 

where  the  superscript  denotes  conjugate  transpose.  Since  the  elements  gPlq(t) 
are  generally  complex- valued,  then  each  G(f)  is  specified  by  2nR  real- valued 
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Figure  3.1:  Structure  of  a  transmitter  with  space-time  encoder  and  time 
varying  linear  transformation. 

parameters  of  which  n|.  are  prescribed  by  (3.3)  and  the  rest  are  free,  but 
constrained  |<?P,9(t)|2  <  ||g9(t)||2  =  1.  These  requirements  can  be  met  by 
elements  of  the  following  form: 

9P,q(t)  -  COS  aPyq(t)e3^o(t\ 

where  aPtq(t),  /3p  q(t)  €  [0, 2ir).  Of  these,  are  independent  realizations 
picked  from  the  uniform  distribution  over  [0, 2ir)  and  the  other  n\  terms  are 
found  from  (3.3).  Of  particular  interest  in  this  report  is  the  case  nT  —  2,  for 
which  an(t),  /3n(t),  /312(t),  /321(t)  are  picked  and  (3.3)  is  solved  for  the  rest 
to  obtain 


G(t)  = 


cos  au(t)e^nW  sin  an^e^1^ 

—  sino:xi(t)e^21^  cosa1i(t)e:?Oi8n(0+/3i2(t)+^2i(t)) 


(3.4) 


3.2  Turbo-STC  with  TVLT 

In  this  section,  properties  of  turbo  space-time  codes  with  TVLT  are  investi¬ 
gated.  It  is  shown  that  full  diversity  can  be  practically  assured.  In  addition, 
if  the  minimum  Hamming  distance  between  all  pairs  of  space-time  codewords 
is  sufficiently  large,  the  coding  gain  exceeds  a  prescribed  threshold  with  high 
probability.  The  Hamming  distance  between  two  codewords  is  defined  as  the 
number  of  columns  that  are  different  between  the  codewords.  Finally,  the 
throughput  rate  of  these  codes  is  investigated  and  is  shown  to  exceed  the 
maximum  rate  for  conventional  space-time  codes  derived  in  [1]. 
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3.2.1  Diversity  and  Coding  Gains 

Let  two  distinct  codewords  of  a  space-time  code  be  represented  by  the 
matrices  C  =  [c(l)...c(L)]  and  B  =  [b(l)...b(L)],  respectively.  According  to 
the  rank  criterion,  the  code  achieves  full  transmit  diversity  if  for  every  pos¬ 
sible  pair  of  codewords  C  and  B,  the  error  matrix  E  =  [e(l), e(L)], where 
e(t)  =  c (t)  —  b (t),  has  rank  rij  [1].  For  special  space-time  codes,  simple 
proofs  of  full  diversity  can  be  found  [1],  however,  for  the  general  case,  as¬ 
certaining  full  diversity  is  a  cumbersome  process  that  requires  testing  the 
rank  of  all  pairs  of  codewords.  When  the  test  fails,  the  encoder  is  redesigned 
and  the  test  starts  over.  For  turbo  space-time  codes,  there  are  no  explicit 
guidelines  for  designing  full  diversity  codes.  According  to  [8],  without  the  te¬ 
dious  check-adjust  process,  a  turbo  space-time  code  with  a  randomly  picked 
interleaver  usually  fails  to  achieve  full  diversity. 

For  a  full  diversity  space-time  code,  the  coding  gain  achieved  by  the  pair 
B,  C  is  £  =  (Ai, ...,  A„T)1//riT  =  (det  (EE^))1^"7  ,  where  Ai,...,A„t  denote  the 
nT  eigenvalues  of  EE*  [1].  The  coding  gain  of  the  space-time  code  is  the 
minimum  of  £  sought  over  all  pairs  of  codewords. 

Consider  now  a  space-time  code  with  TVLT.  Let  two  codewords  be  B  = 
[G(l)b(l), ...,  G(L)b(L)]  and  C  =  (G(l)c(l), ...,  G(L)c(L)] ,  respectively, 
and  the  corresponding  error  matrix  E  =  [G(l)e(l), ...,  G(L)e(L)] .  Let  dh 
be  the  Hamming  distance  between  B  and  C,  i.e.,  the  number  of  nonzero 
columns  in  E.  It  follows  immediately  that  dh  is  also  the  number  of  non-zero 
columns  in  E,  except  for  the  rare  occasion  when  e(t)  =  G(l)e(t)  =  0,  but 
e(t)  7^  0.  Since  the  elements  of  the  transformations  {G(t)}f=1  are  realizations 
of  random  variables,  with  high  probability,  for  dh  >  nT,  rank  E  =  nT,  thus 
giving  the  code  full  diversity. 

The  linear  transformations  {G(t)}f=1  affect  the  diversity  gain  and  coding 
gain  in  a  manner  specified  by  the  following  theorems. 

Theorem  1  Let  B  and  C  be  a  pair  of  nj  x  L  codewords  with  Hamming 
distance  dh ■  Let  d2nax  be  the  largest  Euclidean  norm  of  all  columns  of  E, 
i.e.,  d^ax  =  max  {||e(t)||2}  ,  and  d^in  be  the  smallest  norm  of  all  nonzero 

columns  of  E,  i.e.,  c^in  =  min  {||e(t)||2,  ||e(£)||2  ^  0} .  With  B,  C  the 

£=1,...,L 

corresponding  pair  of  codewords  following  TVLT,  define  the  error  matrices 
E  =  C  —  B  and  E  =  C  —  B.  Then  Pr  jdet  ^EE^  <  u| ,  the  probability  that 

det  (EEf)  is  lower  than  a  given  threshold  u,  where  0  <  u  <  d2nax  •  (dfnin)"r  1 , 


33 


decreases  exponentially  with  [4/  {nT  —  1)J ,  where  [xj  denotes  the  largest 
integer  less  than  x. 

Proof.  Without  loss  of  generality,  rearrange  E  such  that  the  nonzero 
columns  are  listed  first.  In  particular,  let  the  first  column  be  the  one  with 
the  largest  norm  d2nax .  Then,  following  TVLT,  we  have 


E  =  [G(l)e(l),...,G(4)e(4),0,...,0] 
=  [e(l),..,e(4),0,...,0]. 


Divide  the  columns  {e(2),  ...,e(4)}  into  U  —  [ dh/(nT  —  1)J  disjoint  parti¬ 
tions,  each  containing  n^  —  1  columns,  and  concatenate  each  partition  with 
the  first  column  of  E  to  construct  U  square  matrices  V(ju),  p  =  1, ...,  U.  Each 
matrix  V  (p)  is  generated  by  a  set  of  (nT  —  1)  transformations  G  ( t ) .  Since 
the  first  column  is  common  to  all  VQu),  its  transformation  can  be  absorbed 
in  the  MIMO  channel.  Recalling  that  each  G  ( t )  is  obtained  from  n\.  real¬ 
izations  of  a  random  variable,  it  follows  that  its  domain  is  a  hypercube  with 
2 

volume  (27r)"T  .  Extending  this  reasoning  to  V(/i),  its  domain  is  a  hypercube 
with  volume  =  (27r)^nT_1^7lT. 

The  determinant  det  (V (/r)t V (p))  is  a  continuous  function  of  the  angles 
forming  the  transformations  G  (t) ,  since  its  value  is  determined  by  trigono¬ 
metric  and  exponential  functions  of  these  angles.  Moreover,  det  (V (yu)l  V (/i)) 

Tl'f' 

is  bounded  by  0  <  det  (V (/x)fV (p))  <  (V (p)r V {h))kk  -  where  (V (pt)tV (p))  kk 

k= 1 

is  an  element  on  the  diagonal  of  the  matrix  V  (p^  V  (p).  The  left  hand  side 
of  the  inequality  follows  from  the  non-negative  property  of  a  determinant  of 
a  Hermitian  matrix,  and  the  right  hand  side  is  a  consequence  of  Hadamard’s 

Tlx 

inequality.  Direct  evaluation  shows  that  for  any  p,  n  = 

fc=i 


/i(nT-l) 

^max  n  1 1  e  (t )  1 1 2 .  As  the  TVLT  transformations  G  ( t )  range  through 

«=(/i-l)(nT-l)+l 

all  possible  values,  det  (V (p^V (p))  maps  onto  all  values  between  the  two 
bounds. 

[i-(nT- 1) 

Consequently,  for  any  u,  with  0  <  u  <  d?max  {d2min)nT~1  <  d2max  JJ  | |e(t) 1 12, 

t=(/x-l)-(nT-l)+l 

there  exists  a  non-empty  region  in  such  that  if  the  domain  of  V  (p,), 
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A^  £  T^,  then  det  (V(/i)tV(//))  >  u.  Let  denote  the  complement  region 
of  T^,  that  is  if  £  YM,  then  det  (V(p)^V(p))  <  u.  Then,  if  the  volume  of 
TM  is  denoted  |T;(|, 

Pr  {det (V(//)tV(^)  <  u}  (3.5) 

=  ffijl  ...<i 

(27r)("T-1}n| 


Next  the  inequality  det  >  det  (V(/i)^V  (/u))  ,  fi  =  is 

proved  for  later  use.  Since  matrix  E  incorporates  V (p)  and  has  additional 
columns,  denoted  <pl, ...,  4>dh-n-  According  to  the  rank-one  updates  theorem 
for  determinants  [25],  det  (A  +  cd^)  =  det  (A)  (l  +  cP A-1c)  where  A  is 
n  x  n  and  c,  d  are  n  x  1.  Then, 


det  (EEf) 


dh~n 


det  V(//)V(^)t  +  Y  <t>i<t>i 


1=1 


Apply  the  theorem  one  column  at  a  time 


det 


det  ^V(/i)V(/i)f+  Y  <t>i<t>\  j 

det  fv(M)V(M)j  +  Y  &$)  (*  +  0i0i)  ■ 


Repeating  this  step  for  all  columns,  the  claim  is  proved 


det 


=  det 
>  det 


(vWVM'j  n 

i=  1 

(v(  v(rt') . 


(l  + 1 A>!) 


(3.6) 
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With  the  inequality  (3.6)  proved,  it  follows  that, 


Prjdet  (EEf)  <  u} 

<  Pr  {det  (V(l)^V(l))  <  u, ...,  det  (V  (E/)f  V  ([/))  < 

u 

=  npridet  (v(M)tv^))<-} 

ti=i 

=  TT _ EiJ _ 

(27r)(nT“1)nT2 

<  rl^h/(«T-  1)J^ 


where 


r  =  max 
n=i,...,u 


(2^)(nT-1)4  j  • 


(3.8) 


In  (3.7),  line  3  follows  since  (V(/i),/i  =  1,  •  •  •  ,  U}  are  independent;  in  line  5, 
recall  that  U  =  [dh/(nT  —  1)J .  ■ 

The  following  two  theorems  follow  directly  from  Theorem  1. 


Theorem  2  The  probability  that  E  does  not  have  full  rank  decreases  to  zero 
exponentially  with  \dh/{nT  —  1)J  . 


Proof.  According  to  Theorem  1, 

Pr  (rank  E  <  n^} 

=  Pr  | det  {eE*^  =  o} 

<  Pr  | det  (EE*)  <  n}  , 

where  0  <  u  <  d^ax  •  (d^in)nr_1  is  a  threshold  such  that  if  det^EE*^  <  u 
then  the  matrix  EE*  is  claimed  to  lose  full  rank.  Since  Pr  jdet  {EE*^  <  u  j 

decrease  to  zero  exponentially  with  [<4/(nr  —  1)J  ,  then  Pr  {rank  E  <  ut\ 
also  decrease  to  zero  exponentially  with  [dh/{nr  —  1)J  .  ■ 
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Theorem  3  Let  the  coding  gain  achieved  by  the  pair  of  codewords  B,  C 
be  £  =  (det  (EE1))17"  ,  where  E  =  B  —  C.  Then  Pr{£  <  w}  ,  the  prob¬ 
ability  that  the  coding  gain  is  less  than  a  threshold  w,  where  0  <  w  < 

^max  '  (^min)nr  )  decreases  to  zero  exponentially  with  [dh/(nT  —  1)J  . 

Proof. 

Pr  {£  <  w} 

=  Pr{£nT<wnr} 

=  Prjdet^EE1)  <  u}  , 

where  0  <  u  =  wnT  <  (Pmax  ■  (dfnin)nT~1 .  Then  according  to  Theorem  1, 
Pr{£  <  w}  decreases  to  zero  exponentially  with  [dh/(nT  —  1)J  .  ■ 

Theorem  1,  2,  and  3  answer  the  question:  what  is  the  probability  that 
the  code  lose  diversity  and  the  coding  gain  falls  below  a  set  value.  These 
theorems  state  that  both  probability  vanish  exponentially  with  the  ratio 
[ dh/(nT  —  1)J .  It  follows  that  a  desired  feature  of  space-time  codes  with 
TVLT  is  a  large  Hamming  distance.  If  for  all  pairs  of  codewords  the  minimum 
Hamming  distance  is  large  then  the  space-time  code  achieves  full  diversity 
and  good  coding  gain  with  a  high  probability. 

As  to  the  question  how  large  dh,m in  should  be,  there  is  no  general  answer. 
An  example  is  provided  below  to  demonstrate  the  theorems. 

Example  1.  Consider  the  design  of  a  full  diversity  TVLT  space-time  code 
for  nT  =  2  and  4-PSK  modulation.  We  want  to  show  how  large  the  Hamming 
distance  (4,m;n  should  be  to  guarantee  that  the  probability  of  full  diversity  is 
not  achieved  and  the  coding  gain  falling  below  a  specified  value  is  very  small. 
Say  that  Pr  jdet  ^EE^  <  u  j  <  10-6  is  required  for  some  u.  According 

to  (3.7),  if  r  can  be  found,  then  the  required  dh  can  be  determined.  In  this 
example,  with  nT  =  2,  V(/i)  in  (3.5)  has  only  two  columns 

V(p)  =  [e{l),G(p)e(p)\, 

where  without  loss  of  generality  it  was  assumed  that  G  (1)  is  the  identity  ma¬ 
trix  (one  transformation  can  always  be  factored  out  of  the  codeword  and  can 
be  made  part  of  the  channel).  Let  two  pairs  of  transmitted  two-tuple  sym¬ 
bols  with  differences  e(l)  and  e  (/i) ,  respectively.  Also  let  G(/i)  be  given  in 
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(3.4)  and  its  four  angles  taking  values  on  a  grid  in  [0,  2tt)  with  spacings  7r/50. 
For  the  given  symbols  and  for  each  G(p),  we  calculate  det  (V(//)tV(/i))  .  We 
count  the  number  of  realizations  of  G  (p)  that  make  det  (V (p^V (p))  <  u2, 
and  denote  that  number  K  (m,  p,  u ) ,  where  m  indexes  the  pairs  of  two-tuples 
and  p  indexes  the  partition.  Then  the  value  |TM|  in  (3.5)  can  be  approxi¬ 
mated  by  this  K.  We  have 


\r,\  _  K  {m,n,u) 

(27r)(r-l)-»1r  ~  1004 

Checking  over  all  pairs  of  symbols  and  noting  the  largest  count  for  which  the 
coding  gain  fails  the  threshold,  we  have 


K  (p,  u )  =  max  { K  (m,  p,  u)} 

m 


Next,  the  largest  value  of  K  (p,  u)  is  found  over  all  partitions  p  for  the  given 
code.  Let 

K  ( u )  =  max  {K  (//,  u)}  , 

then  from  (3.8),  we  have 


r 


K(u) 
1004  ' 


Assume  a  threshold  of  u  =  d2miJ 8  =  Es/ 4,  where  Es  is  the  energy  of  a 
4-PSK  symbol,  is  set  such  that  if  det^EE^  >u  then  EEt  is  ensured  to  be 
full  rank.  It  is  found  that  r  ~  0.0462.  It  follows  that  a  code  with  a  Hamming 
distance  of  dh,m\n  —  5  makes  Pr  jdet  ^EE^j  <  Es/4  j  «  2  x  10-7  <  10-6. 

In  another  example,  if  u  =  (<^nin)2  /2  =  (Es)2 ,  or  a  coding  gain  threshold 
of  Es  is  required,  it  is  found  that  r  ~  0.27,  and  d^.min  =  11  is  required  for 
Pr  {w  <  Es}  <  10~6.  ■ 

Conventional  space-time  trellis  codes,  such  as  proposed  in  [1],  generally 
do  not  satisfy  a  large  Hamming  distance  requirement.  That  means  that 
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application  of  TVLT  to  such  codes  could  potentially  result  in  low  coding 
gain.  Hence,  TVLT  is  not  suitable  to  be  used  with  conventional  space-time 
trellis  codes. 

Although  not  much  is  known  about  efficient  methods  to  calculate  the  min¬ 
imum  Hamming  distances  for  arbitrary,  non-binary  turbo  space-time  codes, 
such  codes  usually  feature  relatively  large  minimum  Hamming  distances. 
Even  when  there  are  low  weight  codewords  (leading  to  low  Hamming  dis¬ 
tances),  their  number  is  small  and  they  do  not  have  a  significant  impact  on 
performance  [26].  The  calculation  of  for  binary  turbo  codes  is  better 
known.  For  example,  in  [27]  it  is  shown  that  for  a  rate  1/2,  eight-state  binary 
turbo  code  with  convolutional  component  codes  and  interleaver  size  1000,  the 
upper  bound  on  the  minimum  Hamming  distance  is  about  100. 

Hence,  Theorem  1,  2,  3  and  example  1  suggest  that  with  TVLT,  full  di¬ 
versity  as  well  as  good  coding  gain  of  turbo  space-time  codes  can  be  virtually 
ensured. 


3.2.2  Rate 

According  to  Theorem  3.3.1  and  Corollary  3.3.1  in  [1],  for  a  n t  transmit 
and  nR  receive  antennas,  and  signal  constellation  size  26,  the  full  diversity 
of  nTnR  can  be  achieved  with  a  maximum  rate  of  b  bits/s/Hz.  Inspection 
of  the  proof  in  the  reference,  reveals  that  the  rate  limit  is  a  consequence  of 
the  restriction  of  constant  modulation  during  the  codeword.  With  TVLT, 
this  restriction  is  removed,  since  the  modulation  constellation  changes  each 
time  interval.  Consequently,  a  new  rate  derivation  is  required  for  space-time 
codes  with  TVLT. 

The  next  theorem  completes  the  characterization  of  space- time  codes  with 
TVLT,  by  establishing  the  code  rate  and  its  relation  to  the  Hamming  dis¬ 
tance. 

Theorem  4  Consider  a  space-time  code  for  nT  transmit  antennas.  Assume 
a  signal  constellation  Q  (before  TVLT  application)  with  cardinality  |Q|  =  2b 
dements.  Let  L  denote  the  codeword  length.  Then  the  maximum  rate  of  a 
code  with  minimum  Hamming  distance  dh,m in  is 

R  <  nTb  —  -j-  log2  V, 
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in  bits  per  channel  use,  where 


(2 nTb  -  l)J . 


Proof.  Regard  the  columns  of  the  nT  x  L  space-time  codeword  ma¬ 
trix  as  supersymbols  defined  over  a  superalphabet  constellation  Qnj'  with 
\QnT\  —  2”rfc.  Each  codeword  can  be  viewed  as  an  L-tuple  in  a  L-dimensional 
space  [QnT]L .  The  sphere-packing  bound  [28]  can  be  applied  to  obtain  the 
maximum  rate.  Imagine  a  ball  of  radius  m in  in  the  L-dimensional  space, 
centered  at  a  codeword  point,  including  all  points  at  Hamming  distances 
0, 1, ...,  dhtm\n.  The  number  of  such  points  is  given  by  [28] 


v'=  E  ( i 

3=0  '  J  ' 

Associate  one  such  ball  with  each  codeword.  The  total  number  of  L-tuples 
is  (2n'rb)L.  Then  the  maximum  number  of  valid  codewords  separated  by  at 
least  <Li  min  is 

^2  nrb'jL 


A  = 


V 


Hence  the  maximum  transmission  rate  rate  is 

R  <  ~  log2  A 

=  nTb  —  —  log2  V. 


Example  Assume  a  turbo-STC  with  ny  —  2,  L  =  1024,  6  =  3  (8-PSK 
modulation),  and  m;n  =  12.  Then  according  to  Theorem  4,  the  maximum 
rate  is  5.84  bits/channel  use,  which  for  PSK  corresponds  to  the  same  number 
of  bits/sec/Hz.  According  to  [1],  with  constant  modulation,  the  maximum 
rate  for  a  full  diversity  code  is  3  b/s/Hz.  Later  in  this  chapter,  a  rate  4  b/s/Hz 
code  with  8PSK  modulation  that  achieves  full  diversity  is  demonstrated.*  . 

This  section  on  turbo-STC  with  TVLT  is  concluded  with  a  brief  discus¬ 
sion  on  implementation  complexity.  At  the  transmitter,  one  extra  matrix 
multiplication  is  needed  at  every  symbol  interval.  The  added  complexity  is 
not  significant,  especially  when  the  number  of  transmit  antennas  is  two.  Also 
at  the  receiver,  an  extra  matrix  multiplication  is  required  for  every  symbol 
interval  when  compared  to  codes  without  TVLT. 
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3.3  Design  and  Convergence  Analysis 

Conventional  space  time  codes  with  full  diversity  and  high  coding  gain  have  to 
be  either  hand  designed  [1]  or  found  by  exhaustive  search  [10].  This  approach 
is  very  inefficient  for  turbo-STC,  since  two  optimal  component  encoders  do 
not  necessarily  guarantee  a  good  overall  turbo  code,  and  the  search  also  has 
to  include  the  random  interleaver.  In  the  previous  section,  it  has  been  shown 
that  turbo-STC  with  TVLT  has  full  diversity  and  also  a  high  probability  of 
coding  gain.  This  serves  as  a  good  starting  point  in  the  design  of  turbo-STC, 
but  an  analysis  tool  for  these  codes  is  also  needed  to  facilitate  "tweaking" 
the  design.  Such  a  tool  exists  for  binary  [13],  non-binary  [16],  and  coded 
modulation  [24],  in  the  form  of  the  EXIT  chart.  Here,  the  EXIT  chart 
method  is  extended  to  the  turbo  space-time  code  with  TVLT  and  is  used 
to  predict  the  convergence  behavior  of  the  turbo  decoder.  This  extension 
greatly  expedites  the  code  design  process. 

It  will  be  shown  that  with  TVLT,  the  code  performance  depends  on  the  nT 
singular  values  of  the  channel  matrix  H  rather  than  the  channel  matrix  itself. 
In  the  rest  of  the  chapter,  codes  for  two  transmit  antennas  are  exclusively 
considered.  For  this  case,  it  will  be  shown  that  with  =  2,  the  EXIT  chart 
can  be  applied  to  find  the  convergence  region  defined  on  the  2-D  plane  of 
the  channel  singular  values.  Finally,  the  method  by  which  the  chart  can  be 
applied  to  predict  the  code  performance  is  presented. 

3.3.1  Encoder 

Generally  speaking,  space-time  codes  are  multidimensional  codes  [17,29]  with 
the  additional  constraint  of  the  rank  criterion.  With  TVLT,  full  diversity  is 
guaranteed  without  appealing  to  the  rank  criterion.  The  turbo-STC  is  de¬ 
signed  using  accepted  design  practice  for  turbo  codes.  In  Figure  3.2,  is  shown 
the  general  structure  of  a  turbo  space-time  code  encoder  with  TVLT.  The 
component  encoders  are  two  identical,  recursive,  systematic  convolutional 
(RSC)  encoders.  The  output  parity  bits  of  the  two  RSCs  are  alternatively 
punctured.  The  mapping  from  the  coded  bits  to  symbols  is  the  4-D  set 
partition  proposed  in  [29] . 

Several  new  turbo  codes  are  introduced  in  this  report  for  which  the  com¬ 
ponent  encoders  are  listed  in  Table  3.1.  The  notation  used  in  the  table 
and  Figure  3.2  follows  [17].  The  quantities  m  and  v  denote  the  number  of 
input  information  bits  and  the  constraint  length,  respectively.  Each  cir- 
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cle  in  Figure  3.2  represents  a  switch.  The  binary  variable  next  to  each 
switch,  {h\,  i  =  j  =  0,...m}  or  {fk,  k  =  l,...,t>}  ,  indicates  whether 

the  switch  is  in  the  "on"  position.  For  compact  representation,  the  polyno: 
mial  notation  W  (D)  is  used  to  indicate  the  connections  of  the  input  bit  j 
to  the  mod  2  adders  of  the  convolutional  code, 

V 

W(D)  =  J2rtDi. 

i=i 

The  feedback  is  represented  by  H°  (D)  and  the  feedforward  by  F(D), 

F(D)  =  /*(£>). 

fc=l 

All  numbers  shown  in  the  table  are  octal.  For  example,  for  Code  3,  the 
H 3  ( D )  coefficients  are  (04)g  =  (001000)2.  Then 

H3  (D)  =  1  •  D3  +  0  ■  D2  +  0  •  D1  +  0  •  £>°, 

or,  h\  =  0 ,h3  —  0,  hi  =  0,  and  h\  =  1. 

Table  3.1:  List  of  encoders  used  in  this  report,  with  rates  2,  3,  4  b/s/Hz, 
respectively 


Code 

m 

D 

H°(D) 

Hl(D) 

H2(D) 

H3(D) 

H\D ) 

F(D) 

(1) 

2 

m 

68 

18 

28 

38 

(2) 

3 

4 

108 

018 

028 

048 

078 

(3) 

4 

4 

108 

018 

028 

048 

108 

038 

3.3.2  Decoder 

The  structure  of  the  turbo  space-time  decoder  is  shown  in  Figure  3.3.  It  is 
based  on  [7]  and  [11].  It  consists  of  two  MAP  decoders,  which  exchange 


42 


Figure  3.2:  Recursive,  systematic,  convolutional  component  encoder  for 
turbo-STC. 
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extrinsic  information.  The  major  difference  with  the  binary  turbo  decoder 
is  that  the  information  passed  from  one  component  decoder  to  another  may 
be  a  mixture  of  extrinsic  and  systematic  information,  as  shown  next.  In 
binary  turbo  decoders,  the  information  exchanged  between  component  de¬ 
coders  contains  only  extrinsic  information. 

Let  L p  and  La  denote  respectively,  the  a  posteriori  and  a  priori  infor¬ 
mation  in  log-likelihood  form.  Consider  any  one  of  the  component  decoders. 
Then  there  are  two  possible  cases:  (i)  at  time  t  the  parity  bits  from  the  corre¬ 
sponding  component  encoder  are  transmitted;  (ii)  at  time  interval  t  the  parity 
bits  from  the  encoder  are  punctured.  Then  the  a  posteriori  information  can 
be  expressed 


t  (f\  _  (  ka  (t)  +  Lei  (t) ,  (i) 

p[  )~  1  La  (t)  +  Le  (t) ,  (ii) 


(3,9) 


where  Le  (t)  denotes  the  extrinsic  information  and  Le.s  (t)  denote  the  mixture 
of  extrinsic  and  systematic  (intrinsic)  information.  The  expression  for  the  a 
priori  information  is: 


l“W={  liS:©  •  (3-io> 

For  more  details,  see  the  reference  above.  These  differences  with  binary 
turbo  decoder  lead  to  differences  in  the  generation  of  the  EXIT  charts  for 
turbo-STC  as  discussed  in  Section  3.3.4. 

3.3.3  TVLT  Channel 

In  this  subsection,  the  effect  of  TVLT  on  the  channel  model  (3.1)  is  discussed. 
The  singular  value  decomposition  (SVD)  of  the  channel  gain  matrix  H  is 
H  =  USW*,  where  U  is  unitary  nT  x  and  for,  nj-  =  2,  W  is  unitary 
2x2.  The  diagonal  matrix  X  contains  the  singular  values  of  H.  With  that, 
the  channel  model  (3.1)  becomes 


r(f)  =  USWtG(f)c(f)  +  n(f) 

(3.11) 

Ufr  (t)  =  SWfG(f)c(t)  +  U^n(t). 

(3.12) 

Denote  y(t)  =  Ufr(t),  and  Q(t)  =  W^G(i),  with  Q (t)  a  random  unitary 
matrix  that  changes  each  time  interval.  Also,  let  the  noise  term  rj(t)  = 
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U^n(i).  Then  the  equivalent  channel  channel  model  is 

y(t)  =  m (t)c(t)  +  r,(t).  (3.13) 

Since  Q(t)  changes  randomly,  the  equivalent  channel  model  is  a  "forced  fast 
fading"  channel.  When  the  transmitted  codeword  is  infinitely  long,  it  will 
experience  all  possible  values  of  Q  (t).  Consequently,  the  channel  state  is 
parameterized  only  by  S,  or  equivalently  by  the  two  singular  values  oq  and 
a 2.  If  the  codeword  is  long  enough,  the  performance  of  a  code  will  be  invariant 
over  different  channels,  which  have  the  same  singular  values. 

By  applying  TVLT,  the  number  of  parameters  needed  to  characterize  a 
2  x  nR  channel  has  been  reduced  to  one  or  two  real  numbers.  In  contrast, 
without  TVLT,  the  channel  is  described  by  2 N  complex  numbers.  This 
reduction  makes  it  practical  to  apply  the  EXIT  chart  to  the  analysis  of  turbo- 
STC  as  shown  next. 

3.3.4  EXIT  Chart 

The  EXIT  chart  method  [13-16]  is  a  convenient,  low  complexity  tool  used 
to  analyze  the  convergence  of  iterative  decoding  and  to  predict  performance 
of  turbo  codes.  The  turbo  decoder  consists  of  two  component  maximum  a 
posteriori  (MAP)  decoders,  which  exchange  extrinsic  information  as  the  it¬ 
eration  proceeds.  Convergence  to  the  correct  decisions  depends  on  whether 
each  component  decoder  can  generate  better  extrinsic  information  than  it 
receives  from  the  other  decoder.  One  assumption  of  this  method  is  that  the 
codeword  length  or  interleaver  size  are  long  enough  such  that  the  extrin¬ 
sic  information  log-likelihood  ratios  (LLRs)  for  different  time  intervals  are 
independent  [13]. 

The  EXIT  chart  is  a  graphical  representation  of  the  input/output  relation 
of  a  decoder  obtained  through  simulations.  Since  usually  the  two  component 
decoders  are  the  same,  only  one  decoder  needs  to  be  tested.  To  perform  the 
test,  one  needs  to  feed  the  component  decoder  with  the  a  priori  LLRs  of  the 
transmitted  data  symbols  as  if  they  were  passed  from  the  other  decoder.  The 
decoder  is  also  fed  with  channel  observations,  generated  according  to  a  given 
channel  model  H  and  SNR.  The  decoder  is  run  and  the  extrinsic  information 
is  measured  at  its  output.  The  EXIT  chart  is  then  the  plot  of  the  function 

lout  =  f(Iin,  H,  SNR),  (3.14) 
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where  Iin  is  the  mutual  information  between  the  transmitted  data  symbol  and 
the  a  priori  information  (in  LLR  form),  and  Iout  is  the  mutual  information 
between  the  transmitted  data  and  the  output  extrinsic  information  (also  in 
LLR  form). 

The  most  important  application  of  the  EXIT  chart  method  is  to  predict 
the  performance  of  turbo  codes.  For  a  single  antenna  AWGN  channel,  H  can 
be  described  by  a  unit  gain.  Then  the  EXIT  chart  is  the  plot  of  the  function 

lout  =  f  (Iin,  SNR).  (3.15) 

For  a  given  SNR,  tests  for  different  Iin  can  be  done.  And  if  it  is  found  that 

lout  >  hn ,  V/m  >  0,  (3.16) 


then  the  decoder  will  converge  to  correct  decisions  with  high  probability. 
Otherwise,  the  decoding  will  fail.  It  has  been  shown  in  [13]  that  Iout  is  a 
nondecreasing  function  of  SNR.  By  generating  the  EXIT  charts  (3.15)  for 
increasing  SNR  values,  a  threshold  SNRt;,  can  be  found  such  that  if  SNR 
>  SNRift,  then  (3.16)  holds,  and  the  decoder  converges.  The  usefulness  of 
SNRt/j  is  to  predict  the  beginning  of  the  waterfall  region  of  the  bit  error  rate 
(BER). 

For  a  single  input  single  output  (SISO),  block  fading  channel,  the  channel 
7 1  can  be  modeled  by  a  complex  scalar  h  that  follows  some  distribution  and 
changes  independently  between  codewords.  The  EXIT  chart  is  the  plot  of 
the  function 

lout  =  f(Iin,  |h|2SNR). 

The  threshold  SNRt^  found  from  the  EXIT  chart  in  the  AWGN  channel 
can  be  used  to  predict  the  frame  error  rate  (FER)  over  the  block  fading 
channel.  The  convergence  region  for  the  SISO,  block  fading  channel  is 
[h  |  \h\2  >  SNR^/SNR} .  The  the  frame  error  rate  is  the  probability  that 
the  channel  gain  falls  into  the  region  where  the  decoding  fails  to  converge, 


FER  -  Pr 


< 


SNIM 
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of 


For  an  nT  x  nR  MIMO,  block  fading  channel  ,  the  EXIT  chart  is  the  plot 

lout  =  f  (Iin,  HtHSNR). 
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Since  H  is  constituted  from  2  x  x  real  numbers,  the  region  of  con¬ 
vergence  has  to  be  defined  in  a  2  x  rix  x  n^-dimensional  space.  Even  for 
nr  =  2,  Ur  =  1,  searching  the  threshold  boundary  of  a  four-dimensional 
space  is  impractical.  This  fact  impeded  the  development  of  EXIT  charts  for 
turbo  space-time  codes.  As  explained  in  Section  3.3.3,  application  of  TVLT 
solves  this  problem.  The  MIMO  block  fading  channel  is  transformed  into  an 
equivalent  forced,  fast  fading  channel,  which  can  be  parameterized  by  the 
eigenvalues  of  Hi  H  SNR  (or  the  singular  values  of  the  channel  gain  matrix 
H).  For  nr  =  2,  there  are  at  most  two  eigenvalues.  That  is,  the  EXIT  chart 
is  the  plot  of  the  function 

lout  =  /(/in>A  1)A2), 

where  Ali2  are  the  eigenvalues  of  H^H  SNR .  The  EXIT  chart  will  be  used 
to  find  the  region  of  convergence  on  the  2-D  plane  defined  by  the  two  eigen¬ 
values.  The  values  that  lie  outside  the  region  of  convergence  can  be  used  to 
predict  the  FER  performance  of  turbo  space-time  codes  with  TVLT.  Exam¬ 
ples  will  be  shown  in  section  3.3.5. 

To  generate  the  EXIT  chart  for  turbo-STC  with  TVLT,  an  approach  sim¬ 
ilar  to  the  one  proposed  in  Chapter  2  for  SISO  turbo  trellis  coded  modulation 
is  taken.  The  input  to  the  component  MAP  decoder  is  either  generated  as 
pure  extrinsic  information  or  is  synthesized  as  a  mixture  of  extrinsic  and  sys¬ 
tematic  information,  according  to  (3.10).  Also,  the  output  of  the  component 
MAP  decoder  is  either  pure  extrinsic  information  or  mixture  of  extrinsic  and 
systematic  information  according  to  (3.9).  In  the  latter  case,  the  systematic 
information  needs  to  be  removed  before  the  extrinsic  information  is  mea¬ 
sured.  The  synthesis  and  separation  of  extrinsic  information  and  systematic 
information  are  the  major  differences  between  EXIT  charts  for  binary  and 
space-time  turbo  codes.  The  details  of  generating  EXIT  chart  for  multidi¬ 
mensional  turbo  coded  modulation  is  found  in  Chapter  2. 

Figure  3.4  provides  examples  of  EXIT  charts  of  turbo-STC  with  TVLT 
for  three  fixed  and  different  channel  realizations,  but  with  the  same  pair  of 
channel  eigenvalues,  Xi  —  1.5,  A2  =  1.0.  The  code  represented  is  Code  1 
in  Table  3.1.  In  this  figure,  Iin  and  Iout,  measured  in  bits-per-channel-use, 
label  the  abscissa  and  ordinate,  respectively.  The  three  curves  are  plots 
of  I0ut  =  /(/in>  H)  parameterized  by  the  channel  realizations.  Since  each 
transmission  includes  two  information  bits,  /,„  and  Iout  are  in  the  range  [0,2] 
bits.  From  the  figure  it  can  be  seen  that  the  three  EXIT  curves  almost 
overlap  with  each  other. 
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Figure  3.3:  Structure  of  the  turbo  space-time  decoder  with  two  component 
MAP  decoders. 


Figure  3.4:  Demonstration:  turbo-STC  has  same  convergence  behavior  over 
different  channels  with  the  same  eigenvalues. 
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3.3.5  Performance  Analysis 

The  EXIT  chart  is  generated  to  analyze  the  convergence  of  turbo-STC  with 
TVLT  and  to  evaluate  the  FER  performance  over  fading  channel. 


Convergence  Analysis. 

For  the  convergence  analysis  the  following  theorem  is  needed,  based  on 
the  optimality  of  the  constituent  MAP  decoders. 

Theorem  5  Let  Iout  denote  the  mutual  information  between  the  transmitted 
data  symbol  and  the  output  extrinsic  information  from  a  MAP  decoder.  Also 
let  Iin  denote  the  mutual  information  between  the  transmitted  data  symbol 
and  the  input  a  priori  information  to  the  MAP  decoder.  Consider  the  channel 
model  (3.11)  with  nT  =  2  and  £  =  diag(v/Al,  VM),  where  diag(-)  denotes 
a  diagonal  matrix.  Then  I^t  —  f{hn,  Ai,  A2)  is  a  nondecreasing  function  of 
both  Ai  and  A2. 


Proof.  The  channel  model  (3.11)  can  be  regarded  as  an  equivalent  AWGN 
channel  with  symbols  transformed  by  U£WtG(t).  It  will  be  shown  first  that 
in  this  AWGN  channel,  the  Euclidean  distance  between  symbols  is  an  increas¬ 
ing  function  of  the  eigenvalues  Ai,  A2.  Without  loss  of  generality,  consider  a 
scenario  in  which  only  Ai  has  been  increased  to  Ai  +  A,  A  >  0.  Let  d2  (Aj,  A2) 
and  d 2  (Ai  +  A,  A2)  denote  the  respective  Euclidean  distances.  Then  for  V  U, 
W,  ci(i),  c2(t),  Ai,  and  A2, 


d2  (Ai  +  A,  A2)  — 


> 


|| U  diag  (v/ATA,  y/M)  W*G (t)  (d(i)  -  c2(t)) 
||u diag  (0^,  y/%)  W*G (t)  (Cl(t)  -  c2(f))||2 

+  I U diag  (0A,o)  W+G (t)  (c  1(t)  -  c2(t))||2 

||U£1W'fG(t)  (Cl(t)-c2(t))||2 
d2  (Ax,  A2) . 


For  the  scenario  with  the  higher  eigenvalue,  the  MAP  decoder  produces  the 
lower  bit  error  rate  (BER).  If  a  smaller  eigenvalue  was  needed  for  increasing 
the  mutual  information,  then  the  two  scenarios  can  be  made  equivalent  by 
injecting  noise  to  match  the  increase  in  Euclidean  distance.  The  injected 
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Figure  3.5:  The  EXIT  charts  of  the  2  b/s/Hz  code  for  eigenvalues  Ai  =  1 
and  A2  =  0.9,  1.2,  1.5. 

noise  would,  however,  increase  the  BER,  which  contradicts  the  optimality 
of  the  MAP  decoder.  Thus,  MAP  optimality  dictates  that  the  extrinsic 
information  is  non-decreasing  with  the  eigenvalues.  ■ 

Hence,  given  a  Ai,  there  is  a  minimum  A2i  denoted  A2im;n(Ai),  such  that  if 
A2  >  A2imin(Ai),  then  Iout  >  V7j„  >  0,  which  guarantees  the  convergence 
of  the  decoder.  The  value  A2imin(Ai)  is  found  by  gradually  increasing  A2  from 
a  small  value,  until  Iout  climbs  above  the  diagonal  line.  An  example  for  Code 
1  is  shown  in  Figure  3.5.  For  Ai  =  1,  A2,min(l)  is  found  to  be  1.2.  Repeating 
the  process  of  setting  the  value  of  Ai  and  evaluating  A2|min(A1)  by  the  EXIT 
chart,  results  in  a  Ai-A2  characteristic.  In  Figure  3.6,  there  are  shown  Ai-A2 
characteristics  of  the  three  codes  listed  in  Table  3.1.  The  codes  have  rates  of 
2,  3,  and  4  b/s/Hz,  respectively.  For  each  code,  the  area  above  the  curve  is 
the  region  of  convergence.  The  higher  the  data  rate  of  a  code,  the  smaller  the 
convergence  region.  In  addition  to  displaying  the  region  of  convergence,  the 
Ai-A2  characteristic  provides  a  convenient  way  to  compare  and  select  turbo 
space-time  codes  with  TVLT,  by  quickly  predicting  the  FER  performance  of 
the  different  codes,  as  shown  in  the  next  subsection. 
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Figure  3.6:  Convergence  regions  of  the  proposed  codes  with  rate  2,  3,  4 
b/s/Hz,  respectively. 


FER  Analysis. 

The  Ai-A2  characteristics  described  previously  are  applied  to  evaluate 
the  FER  of  turbo-STC  with  TVLT.  For  a  2  x  nR  Rayleigh  fading  channel, 
where  nR>  2,  the  joint  distribution  of  the  eigenvalues  of  H*HSNR  is  given 
by  [30] 


p(Al’ A2)  “  2sSi? exp  (~sfS? (Al  +  A2> (Al  “  A2’2) '  (3'17)- 

Denote  the  region  of  non-convergence  F.  The  FER  of  the  code  at  a  given 
signal  to  noise  ratio  SNR  is  the  integration  of  p(Ai,  A2)  over  the  region  of 
non-convergence  such  as  shown  in  Figure  3.6.  That  is, 


FER  = 


p(Ai,  A2)dA1dA2. 


(3.18) 


The  integration  in  (3.18)  can  be  solved  either  analytically  or  numerically.  The 
accuracy  of  (3.18)  will  be  demonstrated  in  the  next  section  by  comparison 
to  simulation  results. 
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An  analytical  expression  of  the  FER  can  be  found  for  the  2x1  channel. 
In  this  case,  fPH  SNR  has  a  single  eigenvalue.  Let  H  =  [hi  /i2]t,  then 
A  =  (|hi|2  +  |h2|2)  SNR.  For  Rayleigh  fading,  the  eigenvalue  A  follows  a  chi- 
square  distribution  with  four  degrees  of  freedom 


P(x) 


2  SNR2 


(3.19) 


The  minimum  value  of  A  for  which  the  code  converges  is  the  point  (0,  A2)min) 
in  Figure  3.6.  This  value  is  denoted  as  A th,  and  from  the  figure  find  for  Code 
1,  Ajh  ~  3.8.  Then  the  FER  is  expressed  as 


FER  = 


(3.20) 


3.4  Numerical  Results 


In  this  section,  numerical  results  are  provided  for  the  new  codes  with  TVLT 
(Table  3.1)  and  their  performance  is  compared  to  the  case  in  which  no  TVLT 
is  applied.  The  interleaver  length  is  1024  for  all  codes.  The  performance  is 
measured  in  terms  of  the  FER  versus  the  SNR. 

The  performance  of  Code  1,  with  a  spectral  efficiency  of  2  b/s/Hz  with  4- 
PSK  modulation,  is  shown  in  Figure  3.7,  for  both  2x1  and  2x2  configurations 
and  with  and  without  TVLT.  It  is  observed  that  with  TVLT,  the  code  has 
a  steeper  FER  decline  with  SNR,  and  both  2x1  and  2x2  configurations 
are  visually  parallel  with  the  outage  probability  curves.  Without  TVLT, 
the  codes  display  a  marked  loss  in  diversity  gain  in  both  configurations. 
The  turbo-STC  with  TVLT  gap  to  outage  capacity  is  about  1-1.5  dB.  This 
example  demonstrates  the  ability  of  TVLT  to  provide  full  diversity  gains 
(as  discussed  in  earlier  sections).  The  FER  curve  obtained  from  the  Ai-A2 
characteristic  is  also  shown  in  the  figure.  It  is  observed  that  the  prediction 
fits  the  simulation  sufficiently  well  that  can  be  used  as  a  quick  tool  to  evaluate 
performance  of  turbo  space- time  codes. 

The  performance  of  the  8-PSK,  3  b/s/Hz  Code  2  is  presented  in  Figure 
3.8  for  cases  similar  to  those  in  Figure  3.7.  It  is  observed  that  Code  2  loses 
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diversity  gain  in  the  absence  of  TVLT,  but  the  loss  is  not  as  marked  as  in 
the  case  of  Code  1.  A  possible  explanation  is  that  the  8-PSK  modulation 
here  achieves  a  better  rank  of  the  error  matrix  than  the  4-PSK  modulation 
of  Code  1.  In  fact,  this  observation  reinforces  the  argument  for  TVLT,  which 
achieves  the  rank  criterion  through  diversity  of  the  modulation.  Back  to 
the  figure,  it  is  observed  that  the  FER  predicted  from  the  EXIT  chart  is 
reasonably  close  match  with  the  simulation  results. 

The  last  in  the  series  of  performance  analyses  is  that  of  Code  3  displayed 
in  Figure  3.9.  This  code  achieves  a  spectral  efficiency  of  4  b/s/Hz  with  8- 
PSK  modulation.  As  such,  it  breaks  the  rate  limit  of  full  diversity  space-time 
codes  with  fixed  modulation,  which  for  8-PSK  is  3  b/s/Hz  [1].  Conversely, 
codes  at  4  b/s/Hz  cannot  achieve  full  diversity.  This  is  clearly  evident  for 
the  code  without  TVLT  in  Figure  3.9.  With  TVLT,  however,  full  diversity  is 
regained,  as  predicted  by  Theorem  4  in  Section  3.2.  The  EXIT  chart  predicts 
performance  within  1  dB  of  the  simulation  results. 


Figure  3.7:  FER  performance  comparison  of  the  2  b/s/Hz  Code  1  with  and 
without  TVLT. 
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Figure  3.8:  FER  performance  comparison  of  the  3  b/s/Hz  Code  2  with  and 
without  TVLT. 
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Figure  3.9:  FER  performance  comparison  of  the  4  b/s/Hz  Code  3  with  and 
without  TVLT. 
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3.5  Chapter  Summary 

In  this  chapter,  time- varying  linear  transformations  are  applied  to  designing 
full  diversity  turbo  space-time  codes.  TVLT  consists  of  a  unitary  transfor¬ 
mation  which  converts  conventional,  constant  modulations  into  time  varying 
modulations  practically  making  the  rank  criterion  satisfied.  Unlike  space- 
time  and  turbo  space-time  codes  with  fixed  modulations,  which  may  required 
exhaustive  rank  tests  to  ensure  full  diversity,  the  rank  property  is  inherent  to 
TVLT  codes.  Rules  for  the  construction  of  TVLT  have  been  provided  for  an 
arbitrary  number  of  antennas.  It  is  proved  in  a  set  of  theorems  that  the  prob¬ 
ability  that  full  diversity  is  not  achieved  as  well  as  the  probability  that  the 
coding  gain  fails  to  meet  a  prescribed  level  decrease  exponentially  with  the 
minimum  Hamming  distance  of  the  code.  These  theorems  have  two  implica¬ 
tions:  (1)  that  the  design  criterion  for  turbo-STC  is  the  Hamming  distance 
between  codewords  (before  TVLT  application);  (2)  since  turbo-STC  typically 
display  high  minimum  Hamming  distance,  they  are  good  candidates  for  use 
with  TVLT.  A  second  theorem  established  the  rates  achievable  with  TVLT 
space-time  codes.  It  was  shown  that  limits  due  to  constant  modulation  over 
the  codeword  duration  do  not  apply  to  TVLT  and  that  significantly  higher 
rates  are  achievable. 

Beyond  practically  affording  full  diversity  gain,  robust  coding  gains,  and 
high  rates,  TVLT  also  enables  convenient  performance  analysis  and  design 
using  the  EXIT  method.  Previously,  the  EXIT  chart  analysis  has  been  devel¬ 
oped  for  binary  and  multilevel  codes,  but  it  has  been  limited  to  SISO  systems. 
It  has  been  shown  that  in  generating  the  EXIT  chart  for  turbo  codes  with 
TVLT,  the  performance  is  averaged  over  channels  specified  by  the  singular 
values  of  the  channel  matrix.  In  particular,  for  the  two  transmit  antenna 
case,  two  eigenvalues  are  sufficient  to  specify  the  region  of  convergence  of  the 
turbo  code.  It  was  also  shown  that  this  region  of  convergence  can  be  applied 
to  evaluate  the  FER  of  the  code.  Subsequently,  the  EXIT  chart  was  found 
to  be  a  useful  tool  in  the  design  and  analysis  of  turbo  space-time  codes. 

Finally,  various  turbo  space-time  codes  with  TVLT  were  demonstrated 
by  simulations  and  by  FER  analysis.  These  codes  display  all  the  properties 
of  full  diversity,  gain  and  rate  as  predicted  by  the  analysis  and  overall  have 
excellent  performance  with  a  gap  of  only  2  dB  to  outage  capacity. 
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Chapter  4 

Multilevel  Turbo  Space-time 
Coding 


In  the  last  chapter,  the  method  of  using  TVLT  to  design  full  diversity  turbo 
space-time  codes  is  proposed  and  three  codes  for  two  transmit  antennas  with 
rate  2,  3,  4  bits/s/Hz,  respectively,  are  presented.  In  this  chapter,  the  method 
to  design  codes  with  higher  data  rate  without  increasing  the  number  of  trans¬ 
mit  antennas  is  studied.  According  to  [1],  for  a  b  bits/s/Hz  space-time  code 
to  achieve  r  order  transmit  diversity,  the  trellis  complexity  (number  of  states) 
is  at  least  That  is,  the  trellis  complexity  increases  exponentially  with 

the  required  data  rate. 

To  achieve  high  data  rate  with  practical  complexity,  multilevel  codes 
(MLC)  [31]  can  be  used.  With  multilevel  coding,  the  data  bits  are  divided 
into  several  groups  and  each  group  is  encoded  by  an  individual  encoder. 
Consequently,  for  a  given  total  data  rate,  the  data  rate  required  for  each 
individual  encoder  is  lowered.  Then  the  trellis  complexity  of  each  individual 
encoder  is  controlled.  Multilevel  space-time  coding  is  first  studied  in  [1],  and 
later  in  [32]  (only  for  fast  fading  channels). 

To  decode  MLC,  multi-stage  decoding  (MSD)  is  proposed  in  [31].  At  the 
first  decoding  stage,  the  first  level  code  is  decoded  treating  signals  from  higher 
levels  as  interference.  At  the  second  decoding  stage,  the  contribution  of  the 
first  level  is  removed  from  the  channel  observation  while  treating  signals  from 
higher  level  as  interference.  All  the  other  levels  are  then  decoded  in  the  same 
fashion,  sequentially. 

In  this  chapter,  multilevel  turbo-STCs  with  TVLT  are  proposed.  With 
TVLT,  full  diversity  can  be  guaranteed.  Multidimensional  set  partition  tech- 
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niques  are  also  used  in  the  proposed  scheme.  The  advantage  of  multidimen¬ 
sional  set  partition  is  to  increase  the  Euclidean  distance  between  the  modula¬ 
tion  points  seen  by  the  codes  applied  at  higher  levels.  The  scheme  proposed 
in  this  chapter  is  referred  to  as  multilevel  turbo-STC  with  TVLT.  For  simplicT 
ity  of  representation,  in  the  rest  of  this  chapter  the  scheme  is  also  referred  to 
as  multilevel  turbo-STC  without  confusion. 

To  decode  the  multilevel  turbo-STC,  an  iterative  multi-stage  decoding 
(iterative  MSD)  algorithm  is  proposed.  In  conventional  MSD,  only  the  con¬ 
tributions  of  lower  level  codes  are  cancelled.  The  decision  of  the  higher  level 
code  is  not  used  to  help  the  decoding  of  lower  level  code.  In  iterative  MSD 
however,  when  decoding  of  a  certain  level  is  performed  after  the  first  iter¬ 
ation,  the  soft  decisions  of  all  other  levels  from  the  previous  iteration  are 
used  to  remove  their  contributions  from  the  channel  observations.  Conse¬ 
quently  the  system  performance  can  be  improved.  At  the  same  time,  since 
turbo  decoder  is  iterative  anyway,  the  implementation  of  iterative  MSD  with 
turbo-STC  will  not  cause  much  complexity  increase. 

As  an  example  of  multilevel  turbo-STC,  a  two  level,  5  bits/s/Hz  code 
with  16  QAM  modulation  for  two  transmit  antennas  is  proposed.  Simulation 
results  are  also  presented  and  are  compared  to  the  outage  capacity. 

The  rest  of  the  chapter  is  organized  as  follows.  In  section  4.1,  the  encoder 
structure  is  described.  Section  4.2  presents  the  decoding  algorithm.  The 
simulation  results  are  presented  in  section  4.3.  Section  4.4  draws  conclusions. 


4.1  Multilevel  Turbo-STC  Encoder 

The  structure  of  the  multilevel  turbo-STC  encoder  is  shown  in  Figure4.1.  At 
each  time  interval  t,  the  information  bits  are  partitioned  into  W  groups  of 
n* , . . . ,  n'f  bits,  respectively,  and  dispatched  to  W  individual  turbo-STC  en¬ 
coders  to  be  encoded.  In  this  chapter,  the  term  individual  encoder/decoder 
is  used  to  denote  the  turbo-STC  encoder /decoder  for  each  level.  This  term 
should  be  differentiated  from  the  term  component  encoder/decoder  of  which 
each  individual  turbo-STC  encoder/decoder  is  composed.  For  the  multilevel 
turbo-STC  studied  in  this  chapter,  each  individual  turbo-STC  encoder/decoder 
is  composed  of  two  component  encoders/decoders. 

Denote  the  information/systematic  bits  input  to  the  w-th  individual  en¬ 
coder  as  s w(t)  =  [sq  (t),  —,s^_i(t)]  ,  w  =  After  encoding,  individ¬ 
ual  encoder  generates  parity  bits  respectively,  denoted  as 
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Let  uw(t )  =  [s“(t),p“(£)]  denote  the  coded  bits,  including  both  systematic 
bits  and  parity  bits,  from  the  w-th.  individual  encoder.  Further  let  u(t) 

=  ...,  uw (t)  denote  the  coded  bits  from  all  individual  encoders.  The 

coded  bits  v(t)  are  mapped  into  nT  symbols  c(t)  =  [ci(f), ...,  cnT(t)]r  through 
a  given  mapping  M ,  i.e.,  c(t)  —M  (y{t)) .  In  the  proposed  multilevel  turbc^ 
STC  scheme,  M  is  a  multidimensional  set  partition.  With  TVLT  as  de¬ 
scribed  in  Chapter  3,  the  vector  of  symbols  c(t)  are  multiplied  by  a  unitary 
matrix  G(t)  to  obtain  another  vector  x(t)  =  G(t)c(t),  and  are  transmitted 
through  tit  antennas. 

The  major  design  procedure  of  MLC  encoder  is  to  determine  the  rate 
of  each  individual  encoder.  The  rate  design  of  MLC  over  AWGN  channel  is 
studied  in  [31]  [33],  while  the  rate  design  for  MLC  over  the  fast  fading  MIMO 
channel  is  studied  in  [32],  The  basic  idea  is  to  treat  the  MLC  system  as  a 
group  of  individual  codes  transmitting  over  channels  with  different  capacities. 
The  lowest  (first  decoded)  level  has  the  lowest  signal  to  interference  and 
noise  ratio  (SINR)  since  the  interference  from  other  levels  is  not  removed. 
For  higher  levels,  the  interference  from  lower  levels  has  been  removed  before 
decoding  and  the  SINR  is  improved.  Hence  the  capacity  of  the  channel  seen 
by  a  higher  level  is  larger  than  that  seen  by  a  lower  level.  The  code  rates 
for  each  levels  should  then  be  assigned  according  to  these  different  channel 
capacities.  For  the  AWGN  and  fast  fading  channels  considered  in  [31],  [32], 
and  [33],  the  channel  capacities  for  each  level  can  be  easily  computed  since 
the  power  of  interference  from  other  levels  can  be  computed. 

For  the  block  (quasi-static)  fading  MIMO  channel  model  as  studied  in  this 
report,  the  channel  capacity  is  changing  from  one  codeword  transmission  to 
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another.  This  means  the  optimum  rates  for  each  layer  should  also  change 
adaptively  according  to  each  channel  realization.  Adaptive  rate  transmis¬ 
sion  can  be  implemented  in  a  system  with  feedback  channels,  which  is  an 
interesting  future  work  but  is  not  considered  in  this  report. 

In  this  report,  the  rate  design  for  each  level  is  performed  empirically  and 
takes  advantage  of  the  multidimensional  set  partition  (MDSP)  [17].  With 
MDSP,  the  higher  level  has  a  higher  SNR  than  the  lower  level  because  of 
the  enlarged  Euclidean  distance  between  the  modulation  points.  Hence  the 
MIMO  channel  (including  the  modulator)  seen  by  the  higher  levels  have 
higher  capacities  than  the  lower  levels.  Then  the  rate  of  each  level  should  be 
assigned  in  an  increasing  order  from  the  lowest  level  to  the  highest.  But  the 
exact  rate  of  each  level  is  determined  empirically.  Encoders  of  the  higher  level 
can  also  be  designed  with  less  trellis  complexity  while  achieving  comparable 
performance  to  the  lower  level. 

The  method  of  applying  MDSP  to  design  multilevel  turbo-STC  is  demonT 
strated  through  an  example.  A  rate  5  bits/s/Hz,  two  level  (L  =  2)  turbo- 
STC,  with  two  transmit  antennas  ( nT  =  2)  and  16-QAM  modulation  is 
considered.  The  structure  of  the  individual  encoders  is  shown  in  Figure  3.2 
in  Chapter  3.  The  parameters  of  these  two  turbo-STC  encoders  are  shown  in 
Table  4.1.  At  each  time  interval  t,  encoder  1  outputs  systematic  bits  sj(i), 
s\(t),  sl(t)  and  parity  bits  pj(£),  p\(t),  while  encoder  2  outputs  systematic 
bits  sl(t),  sf(£),  and  parity  bits  pg(£). 

The  2-D  set  partition  of  16-QAM  is  shown  in  Figure  4.2.  After  two 
times  partition,  each  subset  Co,  ...,  C3  contains  only  4  points  with  doubled 
minimum  Euclidean  distance  between  points.  Based  on  the  2-D  set  partition 
of  16-QAM,  the  4-D  set  partition  of  2  x  16  QAM  is  shown  in  Figure  4.3.  After 
four  times  partition,  16  subsets  To, ...,  Ti5  are  obtained  with  16  points  in  each 
subset.  It  can  be  shown  that  the  minimum  Euclidean  distance  between  the 
modulation  points  in  each  subset  is  twice  as  large  as  that  of  the  original 
constellation.  In  the  example  code,  four  bits  {Po{t),p\(t),  sj(£),  s](£)}  are 
used  to  select  a  subset  in  {T0, ...,  Ti5}  ,  as  demonstrated  in  Figure  4.3.  The 
remaining  four  bits  {sl{t),Pi{t),  so(0>  si(0}  are  used  to  select  a  point  in  the 
subset.  That  is,  all  the  bits  output  by  encoder  2  are  mapped  into  a  point  in 
one  of  the  subsets  T0, ...,  T15,  which  has  a  twice  as  large  minimum  Euclidean 
distance  as  the  original  constellation  used  by  encoder  1.  Hence,  encoder  2 
should  be  designed  to  have  a  higher  rate,  and/or  less  trellis  complexity  than 
encoder  1.  In  fact,  as  shown  in  Table  4.1,  encoder  1  is  designed  to  be  rate 
3/5  and  each  component  encoder  has  16  states.  In  contrast,  encoder  2  is 
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designed  to  have  rate  2/3  and  each  component  encoder  has  8  states. 


ANT  nT 


ANTI 


Figure  4.1:  Structure  of  the  encoder  of  multilevel  turbo  space-time  code, 
with  set  partition  and  TVLT. 


Table  4.1:  Parameters  of  the  Two  Turbo-STCs  Used  to  form  the  Multilevel 
Turbo-STC 


Code 

Rate 

V 

HU(D) 

H'(D) 

H2(D) 

HA(D) 

H\D) 

F(D) 

Enc  1 

2/3 

3 

68 

18 

28 

0 

38 

Enc  2 

3/5 

4 

10 

01 

02 

04 

0 

1 

4.2  Iterative  Multistage  Turbo-STC  Decoder 

In  this  section,  the  decoding  algorithm  for  multilevel  turbo-STC  is  studied. 
The  channel  model  used  is  the  same  as  (3.1)  and  (3.2)  ,  and  is  reproduced 
below  for  convenience.  The  received  signal  at  time  interval  t  is 

r(f)  =  HG(f)c(t)  +  n(t),  (4.1) 
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Figure  4.2:  Partition  of  16  QAM  constellation  with  increasing  minimum 
subset  distances. 
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Figure  4.3:  Partition  of  a  2  x  16  QAM  constellation  based  on  the  16  QAM 
set  partition. 


where  c (t)  is  the  nT  x  1  vector  of  the  coded  symbols  with  energy  Es/nT\  r (t)  is 
the  nr  x  ur  received  signal  vector;  H  is  the  nr  x  ur  channel  gain  matrix  which 
consists  of  complex- valued  scalars  htj , i  =  j  =  1, ...,  nT,  modeled  as 

zero-mean,  mutually  independent,  identically  distributed  Gaussian  random 
variables  with  unity  variance;  G(f)  is  the  n?  x  i%t  random  unitary  matrix. 
The  noise  term  n(£)  is  modeled  by  zero-mean,  additive  white  Gaussian  noise 
(AWGN)  with  variance  N0/2  per  dimension. 


4.2.1  Iterative  MSD  Decoder 


The  block  diagram  of  the  iterative  MSD  algorithm  is  shown  in  Figure  4.4. 
Corresponding  to  W  individual  encoders,  there  are  W  individual  turbo  de¬ 
coders  at  the  receiver.  Since  the  coded  bits  from  all  individual  encoders  are 
mapped  together  into  a  set  of  symbols  and  transmitted,  the  channel  observa¬ 
tion  must  be  first  decoupled  before  being  passed  to  the  individual  decoders 
to  be  decoded.  The  whole  MSD  decoder  works  in  an  iterative  fashion.  At 
each  iteration,  the  w-th  signal  decoupler  computes  the  decoupled  channel 
observation  for  the  w-th  individual  decoder,  suppressing  the  contribution  of 
the  signals  from  other  levels  by  using  the  a  posteriori  probability  (APP)  of 
the  transmitted  vectors  of  bits  computed  by  their  corresponding  individual 
decoders  in  the  last  iteration.  The  algorithm  for  each  individual  decoder  is 
the  same  as  in  [34]  [7]  with  the  additional  function  of  computing  the  APPs 
for  the  vector  of  coded  bits.  The  rest  of  this  subsection  describes  the  details 
of  the  decoupler  and  decoder  in  one  iteration.  For  simplicity  of  representa¬ 
tion  and  to  avoid  confusion,  the  iteration  index  is  omitted  in  the  following 
expressions. 

Let  h™  =  [60,  ...,6m_i]  denote  a  length  m  binary  vector  that  represents 

m— 1 

a  nonbinary  integer  k,  i.e.,  k  =  bj 2l.  The  decoupled  channel  observation 

i= o 

for  the  vector  of  bits  uw(t)  transmitted  by  the  w-th  individual  encoder  can 
be  represented  as  the  sets  of  log-likelihoods 

{log  (p  (r(f)KM  =  bf ,k  =  0,...,2""+""  -  l}  . 


t  =  1  Let  B%  denote  the  set  of  all  possible  vectors  of  v(t)  such  that 

>'•«  =  bf +"?,  i-e., 


Bt 


{*-(0  I  *'”(«) 


bf+"?}. 
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Figure  4.4:  Structure  of  the  iterative  MSD  for  multilevel  turbo-STC  with  L 
levels. 


Also  define  C™  as  the  set  of  modulated  symbol  vectors  corresponding  to  the 
bits  vectors  in  that  is, 

cr  4  {cm  =x  (>/(«))  i  v(t )  €sn , 


Then 


log  (p  (r(t)K(t)  =  b""+"')) 


(4.2) 


I 


=  log 


=  log 


where 


w 

Y  P  (r(0lc(0  =M  MO))  II  Pa  (MO  =  bfeS+np) 


H  exp 

c(0eCf 


2—1 

i^w 


l 


\ 


V{t)~UN0 'mm 2  +  E  (MO  =  b?+4) 


i=l 

i^w 


La  (i/(t)  =  bf fO)  4  log  (pA  (i/4(t)  =  b^))  , 
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k  =  1, 2r4+np  —  1,1  =  1, W  is  the  a  priori  about  the  bit  vector  transmit¬ 
ted  from  the  i-th  individual  encoder.  In  the  proposed  iterative  MSD  scheme, 
the  a  priori  La  ( 1,1  (t)  =  hks  p )  ,  i  =  1, W,  used  in  the  current  iteration 
is  set  to  be  equal  to  the  APP  generated  in  the  last  iteration.  Let 

LP  (i/(t )  =  b?+”'|  R)  4  log  (p  (y\t)  =  b"!+"i)  |  R)) 

be  the  APP  calculated  by  the  W  individual  decoders  from  the  last  iteration, 
where  the  matrix  R  =  [r(l ),...,  r( N)\  denotes  the  channel  observations  of  the 
whole  codeword.  Then, 

La  (i ,<(t)  =  b”S+”')  =  Lp  (■/(<)  =  b“!+”'|  r)  . 

In  the  first  iteration,  La  (vl(t)  =  b^s+Tlp^  is  set  to  be  zero.  The  method  used 

to  compute  the  APP  of  the  vector  of  coded  bits  vl(t)  is  a  straightforward 
extension  of  the  computation  of  the  systematic  bits  s*(t)  [35]. 

4.2.2  List  Sphere  Decoder 

The  purpose  of  multilevel  codes  is  to  avoid  the  exponentially  increase  in 
trellis  complexity  and  consequently  the  decoder  complexity  when  the  data 
rate  increases.  This  goal  has  been  achieved  by  the  encoder  and  decoder 
algorithm  presented  in  the  last  two  subsections.  However,  the  complexity  of 
the  decoupler,  which  computes  decoupled  channel  observation  for  each  level, 
still  increases  exponentially  with  the  number  of  transmitted  bits  at  each  time 
interval,  as  shown  in  (4.2).  In  (4.2),  for  the  w-th  level,  the  set  Cji’  includes 
all  possible  transmitted  symbol  vectors  c(f),  in  which  the  bits  from  the  w-\h 

Tb^  j  ^ 

level  is  a  specific  combination  bfcs  p.  Then  C™  has  2qw  elements,  where 


w 

<hv  =  (n*  +  np)  • 

ij^w 

That  is,  2qw  times  exp(»)  function  calculation  is  needed  to  carry  out  (4.2). 

To  lower  the  complexity  of  the  decoupler,  the  list  sphere  decoder  (LSD) 
[36]  can  be  used.  The  function  of  the  LSD  is  briefly  introduced  here  while 
the  details  are  in  [36].  Let  C  denote  the  set  of  all  possible  %  x  1  symbol 
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vectors.  Given  the  channel  model  (4.1),  let  c(i)  be  an  unconstrained  ML 
estimate  of  c(f),  i.e., 

c(i)  =  (G(t)tHtHG(t))“1HG(t)r(t), 

the  list  sphere  decoder  search  for  the  set  Clsd (t) 

CLsD(t)  =  |c(f)  |  c (t)  E  arg  min"  (||c(t)  -  c(t)||2}  ,||c(t)  -  c(t)||2  <  82\ 

Ncand  _ 

where  arg  min  { | |c(t)  —  c(t)||2}  denotes  the  set  that  includes  Ncand  nearest 

c  (£)€(? 

symbol  vectors  to  c(£);  8  is  the  search  radius.  The  denotation  Ncanci  is  bor¬ 
rowed  from  [36],  where  it  means  number  of  candidates.  So,  all  the  symbol 
vectors  included  in  the  set  Clsd(1)  are  within  distance  range  52  from  c(f); 
and  there  are  at  most  Ncand  vectors  in  Clsd (£)■  Usually  NCJin(i  is  set  to  be 
much  less  than  2Q,  the  number  of  all  possible  transmitted  symbol  vectors. 
After  Clsd^)  is  found,  the  set  C%  in  the  third  line  of  (4.2)  is  replaced  with 
C™lsd  (t)  i  the  intersection  of  C%  and  Cxs£>(t).That  is, 

CZ,LSD(t)  =  CZnCLsD(t). 

If  CZlsd  ( t )  =  0,  the  decoupled  channel  observation,  or  the  log-likelihood 
in  (4.2)  is  set  to  be  a  default  limit  value.  Since  there  are  fewer  elements 
in  Clsd (t)  than  in  Q.",  the  complexity  of  computing  (4.2)  is  significantly 
lowered. 

The  two  parameters  Ncand  and  8  controls  the  complexity  and  performance 
of  the  signal  decoupler.  If  Ncand  and  8  are  set  to  be  too  small,  few  symbol 
vectors  will  be  included  in  Clsd (i)  to  be  used  in  computing  (4.2),  and  the  de¬ 
coupled  channel  observation  might  be  very  inaccurate.  Thus  the  performance 
of  the  receiver  could  be  significantly  degraded.  Of  course  if  Ncand  and  8  are 
set  to  be  too  large  then  the  complexity  reduction  will  not  be  satisfactory. 
Usually  the  values  of  Ncond  and  8  are  selected  empirically. 


4.3  Simulation  Results 

In  this  subsection  the  simulation  results  for  the  two-level,  5  bits/s/Hz,  16 
QAM  modulation  (before  TYLT)  turbo  STC  for  two  transmit  antennas  are 
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presented.  The  details  of  this  code  have  been  presented  in  section  4.1.  The 
length  of  the  interleaver  used  in  each  individual  turbo  space-time  encoder  is 
N  =  1024.  Two  receive  antennas  are  used  at  the  receiver.  The  channel  is 
block  Rayleigh  fading  as  described  in  (4.1).  The  performance  measured  is 
FER  against  SNR  (dB). 

Both  full  complexity  iterative  MSD  decoder  and  reduced  complexity  de¬ 
coder  with  LSD  are  simulated.  For  the  LSD  used  in  the  simulation,  Ncand 
is  set  to  be  40.  In  contrast,  the  number  of  all  possible  2  x  16  QAM  symbol 
vectors  is  256.  That  is,  the  complexity  of  the  signal  decoupler  is  reduced  to 
about  1/6  of  that  of  the  decoupler  without  LSD.  The  search  radius  S  is  set 
to  be  i/382/SNR  empirically,  such  that  5  decreases  with  SNR.  For  example, 
for  SNR=18  dB,  8  «  2.5. 

The  simulation  results  are  provided  in  Figure  4.5  for  the  proposed  mul¬ 
tilevel  turbo-STC  using  full  complexity  iterative  MSD  decoder  and  reduced 
complexity  decoder  with  LSD,  respectively.  The  outage  probability  curve  is 
also  provided  for  reference.  From  the  figure  it  is  observed  that  using  the  full 
complexity  iterative  MSD  decoder,  full  space  diversity  is  achieved  since  the 
FER  curve  is  parallel  to  the  outage  probability  curve.  And  the  performance 
is  only  2  dB  away  from  the  outage  probability. 

It  is  also  observed  that  the  FER  curve  for  the  reduced  complexity  decoder 
with  LSD  diverges  from  that  of  the  full  diversity  decoder,  showing  that  full 
diversity  is  not  achieved  at  higher  SNR.  However,  at  FER=0.01,  the  perfor¬ 
mance  degradation  is  only  1.5  dB  compared  to  the  full  complexity  decoder, 
which  may  still  be  considered  a  satisfactory  trade  off  between  performance 
and  complexity. 


4.4  Chapter  Summary 

In  this  chapter,  the  multilevel  turbo  STC  scheme  is  proposed  to  avoid  the 
exponentially  increasing  decoding  complexity.  An  iterative  MSD  decoding 
algorithm  is  also  proposed.  Full  diversity  and  close  to  capacity  performance 
is  demonstrated  by  simulations.  A  reduced  complexity  decoder  with  LSD  is 
also  proposed  to  provide  trade  off  between  the  complexity  and  performance. 
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FER 


Figure  4.5:  FER  performances  of  the  two  level  turbo  space-time  codes  over 
2x2  block  Rayleigh  fading  channel. 
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Chapter  5 

Layered  Turbo  Space-time 
Coding 


In  Chapter  3,  the  turbo-STC  with  TVLT  scheme  is  studied  mostly  for  two 
transmit  antenna  and  low  to  moderate  data  rate  cases.  The  reason  is  that 
with  that  scheme  the  decoding  complexity  increase  exponentially  with  the 
number  of  antennas  and/or  data  rate.  Multilevel  turbo-STC  is  proposed  in 
Chapter  4  to  design  codes  with  a  higher  data  rate  without  exponentially 
increasing  decoding  complexity.  But  when  very  high  data  rate  is  required 
and  a  large  number  of  transmit  antennas  is  used,  the  complexity  of  multilevel 
turbo-STC  is  still  very  high.  In  this  chapter  layered  turbo  space-time  coding 
(LTST)  is  presented  to  design  low  complexity  MIMO  systems  with  a  large 
number  of  transmit  and  receive  antennas. 

Layered  MIMO  systems  avoid  the  high  decoding  complexity  by  avoiding 
joint  decoding  of  the  signal  from  all  transmit  antennas.  In  these  systems 
the  signals  from  each  transmit  antenna  or  each  antenna  group,  referred  to 
as  layers,  are  first  separated  and  then  the  separated  signals  are  decoded. 
To  separate  the  signals  belonging  to  a  given  layer  from  other  layers,  the 
signals  from  other  layers  are  treated  as  interference  and  suppressed  with 
either  linear  or  nonlinear  interference  suppression  techniques.  Various  layered 
MIMO  systems  have  been  proposed  recently  [5]  [37]  [38]  [39]  [40]  [41]. 

The  layered  MIMO  systems  can  be  divided  into  two  categories:  Bell  Lab 
Architecture  of  Space-Time  (BLAST)  and  its  variants  [37]  [38]  [39]  [40]  [41], 
etc.,  and  layered  space-time  coding  (LSTC)  [5],  etc.  The  major  difference  be¬ 
tween  these  two  categories  is  the  antennas  included  in  each  layer.  In  BLAST 
systems  each  layer  contains  only  one  antenna,  and  traditional,  well  developed 
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1-D  codes  can  be  used.  The  major  advantage  of  BLAST  systems  is  simple 
design.  In  LSTC,  each  layer  contains  multiple  (usually  two)  antennas  where 
space-time  coding  is  applied.  Compared  to  BLAST  systems,  LSTC  has  fewer 
layers  to  separate  at  the  receiver  and  for  each  layer  there  is  one  less  interferer. 
With  joint  decoding  the  signals  from  all  the  antennas  within  a  layer,  more 
transmit  diversity  is  also  achieved. 

LTST  is  a  special  type  of  LSTC,  with  turbo-STC  applied  to  each  layer. 
Two  LTST  structures  are  proposed  and  studied  in  this  chapter.  One  is  verti¬ 
cal  LTST  (V-LTST),  the  other  one  is  horizontal  LTST  (H-LTST).  The  terms 
vertical  and  horizontal  are  adopted  from  [42]  to  indicate  how  the  coding  is 
applied  to  the  layered  structure. 

In  V-LTST,  only  one  turbo-STC  encoder  is  used  at  the  transmitter.  The 
output  of  the  encoder  is  divided  into  substreams  and  transmitted  from  dif¬ 
ferent  layers.  At  the  receiver,  the  layers  are  separated  using  low  complexity 
parallel  interference  cancellation  (PIC)  techniques.  V-LTST  is  proposed  for 
MIMO  systems  with  low  order  modulation  such  as  BPSK  and  4-PSK.  It  is 
also  suitable  for  low  delay  communication  systems.  In  H-LTST,  an  inde¬ 
pendent  turbo-STC  encoder  with  TVLT  is  used  in  each  layer  and  successive 
interference  cancellation  (SIC)  is  used  at  the  receiver.  The  TVLT  technique 
is  used  to  enable  ordered  SIC  (OSIC).  H-LTST  is  used  for  systems  utiliz¬ 
ing  higher  order  modulation  such  as  8-PSK,  16-QAM,  etc.  The  details  of 
these  two  schemes  are  presented  in  the  rest  of  this  chapter,  after  which  the 
advantages  and  disadvantages  of  two  structures  will  be  made  clear. 

The  rest  of  this  chapter  is  organized  as  follows.  In  Section  5.1  the  V- 
LTST  scheme  is  presented.  The  H-LTST  scheme  is  presented  in  section  5.2. 
Section  5.3  draws  conclusions. 
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5.1  Vertical  Layered  Turbo  Space-time  Cod¬ 
ing 

5.1.1  System  Model 

Consider  a  MIMO  communication  system  with  nT  transmit  and  n#  receive 
antennas  {tit  x  Ur)  operating  over  a  flat,  Rayleigh  block  fading  channel.  The 
output  of  this  MIMO  channel  is  represented  by  nT  x  1  vectors  r  (t)  G  CriR 
that  can  be  expressed 

r(0  —  Hc(t)  +  rf(t),  (5.1) 

where  c(f)  is  nT  x  1  vector  of  transmitted  symbols  whose  energy  is  Es/nr, 
and  the  Ur  x  ur  matrix  H  consists  of  the  channel  coefficients.  The  (i,j) 
element  of  H  represents  the  path  gain  from  transmit  antenna  j  to  receive 
antenna  i.  The  matrix  H  consists  of  complex- valued  scalars  modeled  as 
zero-mean,  mutually  independent,  identically  distributed  Gaussian  random 
variables  with  unit  variance.  The  nT  x  1  vector  ct  contains  the  symbols 
transmitted  at  time  t.  The  noise  term  rj(t)  is  modeled  by  zero-mean,  additive 
white  Gaussian  noise  (AWGN)  with  variance  N0/2  per  dimension. 

The  proposed  V-LTST  structure  is  shown  in  Figure  5.1  with  the  top  part 
showing  the  transmitter.  The  information  data  is  encoded  by  a  symbol-wise, 
parallel  concatenated  space-time  turbo  encoder,  which  will  be  discussed  in 
more  detail  later  in  section  5.1.3.  Briefly,  this  is  a  systematic  turbo  encoder, 
which  generates  a  codeword  matrix  that  consists  of  two  rows,  with  each  row 
transmitted  through  an  antenna.  One  row  of  the  codeword  matrix  contains 
systematic  symbols,  while  the  other  carries  the  corresponding  parity  symbols. 
With  this  arrangement,  the  signals  received  from  the  two  transmit  antenna 
are  treated  as  a  single  entity  and  do  not  need  to  be  separated.  The  codeword 
is  interleaved  pair-wise  so  that  columns  of  the  codeword  matrix  are  preserved. 
This  ensures  that  the  codeword  symbols  associated  with  each  trellis  branch 
(in  this  scheme  one  systematic  symbol  and  its  corresponding  parity  symbol) 
are  transmitted  in  the  same  time  slot.  The  function  of  the  interleaver  is  to 
decorrelate  adjacent  codeword  symbols  and  avoid  possible  burst  errors  after 
interference  cancellation  at  the  receiver.  Following  interleaving,  the  two- 
row  codeword  is  separated  into  several  substreams,  which  are  modulated 
and  transmitted  from  groups  of  two  antennas.  The  substream  of  the  signals 
transmitted  from  each  antenna  group  is  referred  to  as  a  layer. 

The  LTST  receiver  is  shown  in  Figure  5.1-b.  First,  the  signals  received 
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b.  V-LTST  receiver 


Figure  5.1:  Structure  of  the  V-LTST  transceiver  with  a  single  turbo-STC 
encoder  and  two  antennas  per  layer. 
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from  different  layers  are  separated  by  cancelling  the  estimated  signals  from 
other  layers  (regarded  as  interference).  Then,  the  separated  signals  of  each 
layer  are  re-assembled,  de-interleaved  to  the  original  order,  and  passed  to 
the  maximum  a  posteriori  (MAP)  decoders.  Soft  decisions  computed  by 
the  decoders  are  fed  back  to  improve  interference  cancellation,  and  the  two 
decoders  exchange  extrinsic  information  as  in  an  ordinary  turbo  decoder. 
This  process  repeats  for  several  iterations.  The  algorithm  used  at  the  receiver 
will  be  described  and  analyzed  in  the  next  section.  The  data  rate  of  the 
proposed  scheme  is  rnr  log2q  bits/s/Hz,  where  q  is  the  cardinality  of  the 
symbol  constellation  and  r  is  the  rate  of  the  turbo  encoder. 

5.1.2  Encoder  Design 
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Figure  5.2:  Structure  of  the  LTST  encoder  with  2x4PSK  set  partition  to 
optimize  the  Euclidean  distance. 

In  this  subsection  the  design  considerations  of  the  V-LTST  encoder  is  dis¬ 
cussed.  It  was  proved  in  [43]  that  the  design  criterion  for  space-time  codes 
which  work  at  low  SNR  (<  10  dB)  is  to  maximize  the  minimum  Euclidean 
distance  between  codewords:  the  same  design  criterion  as  the  one  for  tradi¬ 
tional  trellis  coded  modulation  (TCM).  This  design  criterion  motivates  the 
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application  of  the  multidimensional  set  partition  techniques  [17]  [29],  origi¬ 
nally  applied  in  TCM  to  increase  the  minimum  Euclidean  distance,  to  the 
design  of  the  turbo-STC  encoder  in  V-LTST.  In  particular,  the  set  parti¬ 
tion  for  multidimensional  phase  modulation  proposed  in  [29]  is  applied  in 
V-LTST.  The  multidimensional  set  partition  guarantees  that  the  distance 
between  branches  departing  and  entering  the  same  state  in  the  trellis  is  at 
least  twice  as  large  as  the  minimum  distance  between  modulation  points. 
The  encoder  structure  is  shown  in  Figure  5.2.  The  component  convolutional 
encoder  is  hand-picked  and  not  guaranteed  to  be  optimal. 

5.1.3  Decoder 

The  turbo  decoder  used  in  V-LTST  is  similar  to  the  log-MAP  decoder  for 
turbo-STC  [34].  It  is  necessary  to  compute  a  posteriori  probabilities  (APP’s) 
for  both  the  systematic  and  parity  symbols  of  a  codeword.  Given  the  obserT 
vations  for  the  entire  codeword,  the  log-MAP  decoder  outputs  a  set  of  APP’s 
for  each  codeword  symbol  c^,  i  =  1, ...,  nT,  t  —  1, ...,  N,  where  N  is  the  code¬ 
word  length  or  frame  length.  Denote  the  APP’s  for  cti  as  Pj  —  P{<h,i  =  $j\ 
observations),  j  =  1,  where  {sj,  j  —  1,  •••,<?}  is  the  modulation  constel¬ 
lation  set.  The  expected  value  of  symbol  Ct,j 

9 

Bm  =  £*pP;  (5.2) 

3=1 

is  used  in  the  co-antenna  interference  cancellation  process  as  explained  below. 
In  the  first  iteration,  when  no  estimate  from  the  decoder  is  available,  {cti, 
t  —  1, ...,  N,  i  =  1, ...,  nr}  are  set  to  zero. 

5.1.4  Interference  Cancellation 

In  V-LTST,  the  signals  from  different  layers  are  separated  by  cancelling  es¬ 
timated  signals  from  other  layers.  Essentially,  the  n^Tn^R  channel  was 
broken  up  into  nT/ 2  parallel  2TnflR  sub-channels.  Suppose  the  desired  sig¬ 
nal  is  supplied  by  the  layer  with  index  i.  For  simple  notation,  the  time  index 
t  is  dropped.  Also,  in  the  following  discussion,  subscript  i  is  used  to  denote 
signals  associated  with  all  layers  other  than  i.  Then  (5.1)  can  be  rewritten 
as 
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r  =  HjCj  +  H-c-  +  77,  (5.3) 

where  H,  is  the  channel  matrix  for  layer  i,  an  nR  x  2  matrix  composed  of  the 
two  columns  of  H  pertaining  to  layer  i ,  and  Hj  is  the  nR  x  (nT  — 2)  matrix  that 
is  made  up  of  the  rest  of  the  columns  of  H;  Cj  is  the  2x1  vector  of  symbols  of 
layer  i,  and  cj  is  the  (nT— 2)x  1  vector  of  symbols  from  other  layers.  Denote  Cj 
an  estimate  of  cq  computed  using  (5.2).  Interference  cancellation  is  achieved 
by  subtracting,  from  the  received  vector  r,  the  estimated  contribution  of  the 
interference  Hjc q 

y  i  =  r-H-c- 

=  HjCj  +  Hj(cj  —  cqj  +  rj 
=  H^  +  H^  +  77,  (5.4) 

where  y *  is  the  signal  from  layer  i  cleaned  of  interference,  and  eq  =  (cq  —  cq) 
is  decision  feedback  error.  Prom  (5.4),  it  can  be  seen  that  each  layer  is  a 
2  x  nR  MIMO  channel  with  interference.  Let 

Ri  A  E[(H;e;  +  »j)(Hjej  +  tj)1] 

=  (.E.HjHJ  +  AToI.  (5.5) 

be  the  correlation  matrix  of  the  residual  interference  plus  noise,  where  the 
superscript  denotes  complex  conjugate  transpose.  The  relative  strength  of 
the  interference  term  is  given  by  the  normalized  mean  square  error  of  the 
decisions  feedback  p  =  E(\e~hl\2)/ Es,  where  e-ih  l  =  1,  —  2,  is  an  element 

of  the  vector  ej.  Then  the  optimum  branch  metric  used  in  the  decoder  is 

7(Yi.Xi)  =  (y i  -  Hixi)tRf1(yi  -  HiXf),  (5.6) 

where  Xj  are  the  symbols  to  be  tested. 

The  disadvantage  of  using  (5.6)  as  a  branch  metric  lies  wherein  the  high 
complexity  of  inverting  the  matrix  R,,  whose  dimensions  are  nR  x  nR.  V- 
LTST  takes  a  much  simpler  approach  of  spatial  matched  filtering  (SMF). 
Simulation  results  in  section  IV  demonstrate  that  in  spite  of  the  subopti- 
mal  spatial  interference  suppression,  V-LTST  achieves  performance  close  to 
capacity.  Following  SMF,  the  desired  signal  for  layer  i  is: 

y  i  =  Hjy; 

=  H-HiCj  +  H\Hjeq  +  H-?7 
=  A,Cj  +  Bjej  +  H-?7 

=  AjCj  +  Vi  +  Zi,  (5.7) 
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where  A;  =  HjH;  is  a  2  x  2  matrix;  B;  =  H-Hj  is  a  2  x  (n^  —  2)  matrix;  v*  = 
B,ej  and  z;  =  H\r)  are  respectively,  2x1  vectors  of  the  residual  interference 
and  noise  term  post  SMF.  It  is  noted  that  (5.7)  is  now  an  equivalent  2T2R 
channel  model. 

Assuming  residual  interference  plus  noise,  v*  +  z*,  to  be  AWGN,  the 
following  sub-optimum  branch  metric  is  used  in  V-LTST 

7(yi,Xi)  =  (A')-1  (yi  -  -  HjXj),  (5.8) 

where  Aq  is  the  variance  of  the  elements  of  the  vector  v*  -I-  Zj.The  variance  of 
interference  plus  noise  Nq  is  needed  in  (5.8)  for  incorporation  of  the  a  priori 
information  in  the  MAP  decoding,  and  is  computed  as 


Aq  =  E[(vi  +  Zi)(vi  +  Zi)]] 

=  \  ■  tr(BiBt)  +  A0  •  tr(Ai)^  , 


(5.9) 


where  tr(-)  is  the  trace  operation.  And  the  signal  to  interference  noise  ratio 
(SINR)  after  SMF  is 


In  the  first  iteration  no  soft  decision  is  available  to  cancel  the  signals 
from  other  layers  and  p  =  1.  The  initial  SINR,  with  p  =  1  will  be  at  most 
1/  ^2tr(B,Bt)j  ,  which  will  not  decrease  with  the  SNR=  Es/Nq.  If  the  data 
rate  is  high  and  the  modulation  size  is  large,  this  initial  SINR;  can  be  too 
small  for  the  decoder  to  generate  good  enough  soft  decisions,  to  be  usable 
for  cancellation  in  the  next  iteration.  In  this  report  it  is  found  that  V-LTST 
is  only  suitable  for  small  size  modulation  such  as  4-PSK  and  BPSK. 


5.1.5  Space  Interleaving 

In  V-LTST,  The  dominant  error  events  for  layered  systems  are  caused  by  the 
channel  gain  matrix  for  one  layer  being  nearly  "parallel"  (opposite  to  orthog¬ 
onal)  to  that  of  other  layer(s).  When  the  channel  gain  matrix  of  one  layer 
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is  almost  "parallel"  to  that  of  the  other  layer(s),  the  interference  between 
layers  will  be  strong.  Consider  the  SMF  used  in  (5.7).  When  the  channel 
gain  matrix  of  Hj  for  layer  i  is  orthogonal  to  Hj,  or  IIHjHdl  2  =  0,  there  will 
be  no  interference  from  other  layers  at  all.  The  notation  ||  ■  ||2  represents 
the  Frobenius  norm.  Conversely,  if  for  some  channel  realizations  ||H|Hj||2  is 
large,  the  interference  from  other  layers  will  be  strong.  To  improve  the  sys¬ 
tem  performance  in  these  cases,  the  combination  of  antennas  in  each  layer 
during  the  transmission  can  be  alternatively  changed,  which  is  referred  to  as 
space  interleaving.  A  4T4R  system  will  be  used  as  an  example.  Let  H  =  [hi 
h2  1*3  h4]  where  i  =  1,  ...,4  are  columns  of  H.  The  combination  of  anten¬ 
nas  in  each  layer,  as  it  evolves  with  time,  is  illustrated  in  Table  5.1.  As  a 
result  of  space  interleaving,  if  at  some  time  interval  ||HJHj||2  is  large,  it  is 
possible  that  at  another  time  interval  with  another  combination,  ||H]Hj||2  is 
small. 

Table  5.1:  An  Example  of  Space  Interleaving  of  the  Transmit  Antennas  in 
Two  Layers 


time 

layer  1 

layer  2 

Hi 

h2 

1 

ant  1,2 

ant  3,4 

hih2 

h3h4 

2 

1,3 

2,4 

hih3 

h2h4 

3 

1,4 

2,3 

hih4 

h2h3 

4 

1,2 

3,4 

hxh2 

h3h4 

.  .  . 

.  .  . 

.  *  . 

.  .  . 

.  .  . 

Similar  techniques  are  proposed  for  single  antenna  per  layer  systems  such 
as  [37]  [39]  for  different  purposes.  In  those  schemes,  a  signal  from  each  layer 
is  cycling  through  all  the  transmit  antennas.  In  this  fashion,  each  layer  will 
achieve  more  diversity. 

5.1.6  Complexity  Analysis 

In  this  subsection  and  next  subsection  the  V-LTST  scheme  is  going  to  be 
compared  to  existing  schemes  with  the  same  number  of  antennas,  data  rate, 
and  comparable  performance.  The  turbo-greedy  algorithm  [41]  is  a  turbo 
coded  vertical-BLAST  (V-BLAST)  systems,  with  simplified  joint  maximum 
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likelihood  (ML)  demodulation,  referred  to  as  turbo-greedy  algorithm,  of 
the  bits  transmitted  from  all  transmit  antennas.  The  threaded  space-time 
(TST)  scheme  [40]  uses  a  transmitter  structure  similar  to  diagonal  BLAST 
(D-BLAST),  with  a  minimum  mean  square  error  (MMSE)  demodulator  to 
separate  the  signals  from  each  transmit  antennas. 

Now  the  complexity  of  V-LTST  with  that  of  other  MIMO  schemes  with 
similar  performance  is  compared.  The  complexity  of  interference  subtraction 
and  SMF  is  only  0(nT).  Typically,  the  complexity  of  the  V-LTST  receiver 
is  dominated  by  the  turbo  decoder  whose  log-MAP  decoding  has  a  com¬ 
plexity  of  0(u 2log<J),  where  v  is  the  number  of  states  of  the  convolutional 
encoder  and  q  is  the  modulation  size.  In  contrast,  the  complexity  of  ML 
spatial  demodulation  is  0(exp(nr))  and  that  of  its  simplified  version,  the 
turbo-greedy  algorithm,  is  less  than  0(exp(nT)),  but  higher  than  0(n £)  [41]. 
As  another  point  of  comparison,  the  complexity  of  the  threaded  space-time 
(TST)  scheme  [40]  with  MMSE  interference  suppression  is  0{n\).  Thus  for 
nT  =  4,  q  =  4,  u  —  8,  the  complexity  of  V-LTST  is  0(32)  compared  to 
0(256)  for  turbo-greedy  and  MMSE  interference  suppression. 

5.1.7  Numerical  Results 

Numerical  results  for  the  V-LTST  scheme  are  presented  for  4  x  4  and  for  8  x  8 
configurations.  For  4-PSK  modulation,  the  data  rate  for  the  4x4  system  is 
4  bits/s/Hz,  and  for  the  8x8  system  it  is  8  bits/s/Hz.  Both  systems  are 
full  rate.  The  channel  is  assumed  to  be  constant  over  blocks  of  130  channel 
uses.  Consequently,  the  frame  length  for  4  x  4  system  is  260  symbols  (130 
channel  uses  x  2  layers),  and  for  the  8x8  system  is  520  symbols.  The  SNR 
is  defined  Es/Nq,  i.e.,  the  ratio  of  the  total  transmitted  power  to  the  noise 
power. 

To  determine  the  number  of  iterations,  a  simple  stop  criterion  is  used  in 
V-LTST,  where  the  iterations  stop  whenever  two  successive  iterations  yield 
the  same  hard  decision.  For  all  the  V-LTST  simulation  results  presented 
below,  the  average  number  of  iterations  varied  from  3  to  5.  Numerical  results 
shown  for  TST  and  turbo-greedy  algorithms  are  lifted  from  [40]  and  [41], 
respectively. 

Figure  5.3  shows  the  frame  error  rate  (FER)  for  4  x  4  systems.  It  can  be 
observed  that  V-LTST,  TST,  and  turbo-greedy  curves  are  approximately  par¬ 
allel  indicating  the  same  space-diversity  gain.  V-LTST  outperforms  turbo- 
greedy  coding  by  0.5  dB  and  outperforms  TST  by  approximately  1.3  dB. 
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Figure  5.3:  FER  comparison  of  V-LTST  and  other  schemes,  4  x  4,  4  b/s/Hz. 

It  can  also  be  observed  that  the  performances  of  both  V-LTST  and  turbo- 
greedy  are  within  2  dB  range  of  the  outage  probability.  Figure  5.4  shows 
the  FER  comparison  for  8  x  8  systems.  It  can  be  observed  that  V-LTST 
and  turbo-greedy  coding  have  similar  performances  and  outperform  TST  by 
about  1  dB. 


5.2  Horizontal  Layered  Turbo  Space-time  Cod¬ 
ing 

V-LTST  achieves  close  to  outage  probability  performance  with  simple,  low 
complexity  interference  subtraction  to  separate  the  signals  from  different  lay¬ 
ers.  But  V-LTST  is  only  suitable  for  low  order  modulation  systems.  In  this 
section,  H-LTST  is  then  proposed  for  systems  that  require  a  higher  rate  and 
higher  order  modulations. 
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Figure  5.4:  FER  comparison  of  V-LTST  and  other  schemes,  8x8,  8b/s/Hz. 

5.2.1  System  Model 

The  structures  of  the  H-LTST  transmitter  and  receiver  are  shown  in  Figure 
5.5.  At  the  transmitter,  the  nx  transmit  antennas  are  partitioned  into  n^/2 
layers  with  two  antennas  in  each  layer.  An  independent  turbo-STC  encoder 
is  used  in  each  layer.  The  nr/ 2  encoders  may  be  different.  In  this  chapter 
however,  only  the  case  where  the  encoders  for  every  layer  are  the  same  is 
considered.  The  TVLT  technique  is  also  applied  to  each  layer.  That  is, 
the  coded  symbols  are  multiplied  by  a  random  unitary  matrix  before  being 
transmitted. 

At  the  H-LTST  receiver  SIC  is  used  to  separate  the  layers.  That  is,  the 
interference  from  the  previously  detected/decoded  layers  are  cancelled  while 
the  interference  from  other  undetected  layers  are  nulled  by  ZF  filtering.  After 
the  cancel-null  operation,  the  signals  for  the  desired  layer  is  obtained  and 
decoded. 

For  SIC  detector,  the  detection  order  is  important  since  the  layer  with  the 
least  decoding  error  probability  should  be  detected  first  [38] ,  to  minimize  the 
effect  of  error  propagation  and  then  to  minimize  the  overall  error  rate  of  the 
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Figure  5.5:  Structure  of  H-LTST  with  independent  turbo-STC  en¬ 
coder/decoder  for  each  layer. 

whole  system.  In  H-LTST,  TVLT  is  applied  to  each  layer.  Hence  whether  or 
not  the  decoder  of  a  certain  layer  will  converge  to  the  correct  decision  (error 
free)  can  be  predicted  with  high  accuracy,  using  the  method  presented  in 
Chapter  3.  Then  the  layer  that  is  predicted  to  converge  should  be  detected 
first  and  the  optimum  detection  order  can  be  determined. 

In  the  previous  proposed  layered  space-time  coding  (LSTC)  scheme  [5], 
however,  ordered  SIC,  or  OSIC,  is  not  applicable  and  the  layers  are  detected 
in  an  arbitrary  but  fixed  order  independent  of  the  current  channel  condition. 
To  mitigate  the  performance  loss  due  to  lack  of  optimum  detection  order, 
non-constant  power  transmission  is  used.  In  [5],  to  take  advantage  of  the  fact 
that  later  detected  layers  have  a  higher  diversity  order,  the  earlier  detected 
layers  are  assigned  a  higher  transmit  power  than  the  later  layers.  This  non¬ 
constant  power  transmission,  however,  leads  to  an  increased  requirement  for 
the  amplifier  dynamic  range  and  then  the  amplifier  cost. 

5.2.2  Ordered  SIC 

The  channel  model  is  the  same  as  in  (5.1)  and  (5.3),  except  that  the  trans¬ 
mitted  signal  of  each  layer  is  multiplied  by  the  random  unitary  matrix  G (t), 
referred  to  as  TVLT  in  Chapter  3.  Let  c(f)  be  the  nr  x  1  vector  of  coded 
symbols  to  be  transmitted  at  time  interval  t,  each  having  energy  Es/nx- 
Also  let  c#;(f)  be  the  2x1  vector  of  coded  symbols  from  the  k-th  layer. 
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That  is  c (t)  =  cf,...,c ,  where  (■)T  denotes  matrix  transpose.  Let 

T  lr 

Xfc(t)  =  G(t)ck(t)  be  signal  after  TVLT  and  x  =  I  xf, ...,x^t/,2  be  the 
nT  x  1  vector  of  transmitted  signal.  To  simplify  the  representation,  the  time 
index  is  dropped.  The  channel  model  is  then, 

r  =  Hx  +  n, 

where  r  is  the  n^xl  vector  of  received  signal;  H  is  the  Ur  x  nj-  channel  gain 
matrix;  n  is  AWGN. 

Before  the  OSIC  algorithm  is  presented,  several  variables  are  defined.  Let 
H  [k]  denote  the  channel  gain  matrix  of  the  k- th  layer,  which  is  composed  of 
the  2/c-th  and  the  (2k  +  l)-th  column  of  H.  Let  H  [m,  n]  denote  the  channel 
gain  matrix  pertaining  to  the  layers  from  mto  n.  Also  let  H  [m,  n  \  i ]  denote 
the  channel  gain  matrix  pertaining  to  the  layers  from  m  to  n,  except  layer  z. 

Without  loss  of  generality,  it  is  assumed  that  at  the  k- th  stage  of  OSIC, 
layers  1  ,...,fc  —  1  have  been  detected  without  error.  Denote  the  detected 
and  decoded  symbols  as  c  [1,  A  —  1]  =  [cf , ...,  c^_j]  and  the  correspond¬ 
ing  the  signal  after  TVLT  x  [1,  /c  —  1]  =  [xx,  ...,Xfc-i]T .  Then  the  channel 
observations  pertaining  to  layers  k  to  nrj 2  are  obtained  by  removing  the 
contributions  of  layers  1  to  k  —  1: 

r[k,riT/2]  =  r  —  H  [1,  k  —  1]  x  [1,  k  —  1] 

=  H  [k,nT/2]  x  [k,  tit/2]  +  n, 

where  x  [k,  nr/2]  is  the  signal  transmitted  from  the  layers  k  to  nr/2. 

To  detect  the  next  layer,  the  signal  for  the  other  undetected  layers  should 
be  nulled.  Assume  the  i-ih  layer,  k  <i  <nT/ 2  is  going  to  be  detected.  After 
nulling  the  signals  from  other  layers,  the  signals  of  the  i- th  layer  is 

y[i]  =  $Wrk 

=  ^[i}^H[k,nT/2]x[k,nT/2]  +  <&Jn 
=  [z]  x  [i]  +  3>jn 

=  H  [z]  x  [z]  +  z 

where  <&[z]  is  the  projection  to  the  null  space  of  H  [k,  nT/2  \  z]  whose  columns 
are  normalized;  H  [z]  =  [z]  is  the  post  projection  channel  gain  matrix  for 

the  z-th  layer;  z  =  $[z]^n  is  the  transformed  noise  term.  Since  the  columns 
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of  <&[z]  is  normalized,  the  variance  of  the  noise  is  not  changed  after  the 
projection.  The  matrix  3>[z]  is  obtained  by  first  computing 

#[z]  =  I  -  H  [k, nr/ 2  |  z]  (h  [Jfc,  nT/ 2  |  i}1  H  [k,  nT/ 2  |  z])  ^  H  [k, nT/ 2  |  z]1 
and  then  normalize  its  columns. 

To  minimize  the  error  rate,  the  detection  order  should  be  determined. 
That  is,  the  layer  with  the  least  error  probability  should  be  detected  first. 
For  conventional  layered  space-time  codes  [5],  ordering  is  not  applicable  since 
it  is  difficult  to  predict  which  layer,  with  the  post  projection  channel  H  [z] , 
i  =  k, ...,  nr/2,  will  have  the  least  probability  of  error. 

For  turbo  STC  with  TVLT,  whether  or  not  the  decoder  for  a  specific  layer 
will  converge  to  the  correct  decision,  can  be  predicted  with  high  accuracy 
by  checking  whether  the  eigenvalues  (Ai[z],  A2 [i] )  of  [z]'  H  [z]  fall  in  the 
region  of  convergence,  as  shown  in  Figure  3.6.  Then  the  layer  that  is  predicted 
to  converge  should  be  detected  first.  If  for  multiple  layers  the  eigenvalues 
(Ai[z],A2[z])  fall  in  the  convergence  region,  the  layer  that  should  be  detected 
first  is  determined  as  follows.  Given  an  arbitrary  Ai[z],  let  A2,min  (Aj.[z])  be 
the  minimum  required  A2  for  the  z-th  decoder  to  converge.  Then  the  A:-th 
layer  with 

k  =  argmax  (A2[z]  -  A2,min  (Ai[z])) . 

I 

should  be  detected  first.  That  is,  the  layer  whose  eigenvalues  are  the  farthest 
away  from  the  border  between  the  region  of  convergence  and  region  of  non¬ 
convergence  should  be  detected  first. 

5.2.3  Numerical  Results 

In  this  section  the  simulation  results  for  the  proposed  H-LTST  scheme  are 
presented.  In  H-LTST,  TVLT  is  employed  at  the  transmitter  which  in  turn 
enables  the  OSIC.  For  comparison,  the  performances  of  two  other  systems 
are  also  shown.  One  system  does  not  employ  TVLT  and  consequently  can 
only  perform  SIC  without  ordering.  This  system  has  a  similar  structure 
with  the  conventional  LSTC  [5].  Another  system  employs  TVLT  but  still 
performs  SIC  without  ordering  although  ordering  is  available.  By  comparing 
the  performance  of  H-LTST  with  those  of  the  two  other  systems,  the  effects 
of  TVLT  and  OSIC  will  be  clearly  shown. 
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The  turbo-STC  encoder  used  in  all  systems  are  Code  1  from  Table  3.1 
with  interleaver  length  1024.  This  is  a  rate  3/6  encoder  with  8-PSK  modula¬ 
tion.  Simulations  results  will  be  shown  for  four  transmit  antenna  four  receive 
antenna  (4  x  4)  and  eight  transmit  antenna  eight  receive  antenna  (8  x  8)  con¬ 
figurations.  The  data  rates  of  these  two  configurations  are  6  b/s/Hz  and  12 
b/s/Hz,  respectively. 

The  FER  performances  of  the  4x4  systems  are  shown  in  Figure  5.6.  It  is 
observed  that  H-LTST  achieves  more  diversity  than  the  two  other  systems. 
At  FER— 0.01,  H-LTST  outperforms  the  systems  with  no  TVLT  and  no  or¬ 
dering  by  2.5  dB,  and  outperforms  the  system  with  TVLT  but  no  ordering 
by  2.0  dB.  This  comparison  shows  the  importance  of  ordering  in  SIC  detec¬ 
tion,  and  demonstrates  the  advantage  of  H-LTST  over  conventional  LSTC  [5] 
which  employs  SIC  detection  but  ordering  is  not  available.  It  can  also  be  ob¬ 
served  that  the  system  with  TVLT  but  no  ordering  achieves  more  diversity 
than  the  system  with  neither  TVLT  nor  ordering,  but  the  FER  performance 
improvement  caused  by  the  diversity  improvement  is  limited.  The  reason 
is  that  although  TVLT  can  guarantee  full  transmit  diversity  of  each  layer, 
with  four  receive  antennas  in  the  systems  which  already  leads  to  a  high  order 
receive  diversity,  whether  or  not  full  transmit  diversity  is  achieved  will  not 
cause  significant  performance  difference. 

The  FER  performance  of  the  8x8  systems  are  shown  in  Figure  5.7.  In 
this  case,  H-LTST  still  achieves  more  diversity  than  the  other  systems.  At 
FER=0.01,  it  outperforms  the  system  with  no  TVLT  and  no  ordering  by 
4  dB,  and  outperforms  the  systems  with  TVLT  but  no  ordering  by  3  dB-. 
The  more  performance  improvement  in  8  x  8  case  than  in  4  x  4  case  shows 
that  as  the  number  of  transmit  antenna  ,or  layers,  increases,  ordering  in  SIC 
plays  a  more  important  role.  It  can  be  also  observed  that  only  applying 
TVLT  to  each  layer  without  using  OSIC  will  lead  only  limited  performance 
improvement. 


5.3  Chapter  Summary 

In  this  chapter,  two  types  of  layered  turbo  space-time  coding  schemes,  V- 
LTST  and  H-LTST  are  proposed  for  MIMO  systems  with  large  number  of 
transmit  antennas.  The  performance  of  the  proposed  schemes  have  been 
compared  to  existing  layered  MIMO  schemes. 

V-LTST  is  proposed  for  systems  with  low  order  modulation  and  low  data 
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Figure  5.6:  Performance  of  H-LTST,  4  x  4,  6  b/s/Hz,  with  comparison  to 
other  schemes. 
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Figure  5.7:  Performance  of  H-LTST,  8  x  8,  12  b/s/Hz,  with  comparison  to 
other  schemes. 


rate  per  layer.  With  two  antennas  per  layer  structure,  the  number  of  in¬ 
terferes  for  each  layer  is  reduced  compared  to  the  BLAST  systems.  Space 
interleaving  is  applied  to  further  reduce  the  interference  between  layers.  To 
maximize  the  coding  gain,  multidimensional  set  partition  is  applied  to  the 
encoder  design.  With  these  design  techniques,  the  resultant  V-LTST  scheme 
avoids  high  complexity  signal  separation  algorithms  such  as  MMSE  and  ML, 
by  performing  only  low  complexity  PIC  in  conjunction  with  turbo  coding.  It 
is  shown  that  V-LTST  achieves  performance  within  2  dB  of  outage  capacity, 
but  with  significantly  less  complexity  than  other  MIMO  schemes  with  similar 
performances. 

H-LTST  is  proposed  for  systems  with  high  order  modulation  and  high 
data  rate  per  layer.  At  the  transmitter,  TVLT  is  applied  at  each  layer  with 
turbo-STC.  Because  of  the  TVLT  used  at  each  layer,  whether  or  not  the 
decoder  of  each  layer  will  converge  can  be  predicted  based  on  the  eigenvalues 
of  the  channel  gain  matrix  of  each  layer.  Hence  the  detection/decoding  order 
can  be  determined.  Simulation  results  show  that  H-LTST  achieves  significant 
performance  improvement  upon  layered  space-time  coding  without  ordered 
detection /decoding. 

Performance  improvements  by  the  proposed  schemes  have  been  demon¬ 
strated  by  comparing  them  to  other  existing  schemes.  It  is  also  worth  com¬ 
paring  V-LTST  with  H-LTST.  The  major  advantage  of  V-LTST  is  that  it  is 
more  suitable  for  communication  systems  with  short  delays.  It  is  known  that 
turbo-STC  works  well  only  when  the  codeword/interleaver  is  long  enough. 
For  a  system  with  K  channel  use  per  frame,  the  length  of  the  codeword 
transmitted  using  V-LTST  is  Kut/ 2  since  signals  transmitted  from  all  lay¬ 
ers  belong  to  a  single  turbo-STC  encoder.  With  H-LTST,  each  layer  uses  an 
independent  encoder,  the  length  of  the  codeword  transmitted  in  each  frame 
equal  to  K.  The  advantage  of  H-LTST  is  that  it  can  be  applied  to  systems 
with  high  order  modulation  while  V-LTST  is  only  suitable  for  low  order 
modulation  systems. 
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Chapter  6 

Summary  and  Future  Work 


In  wireless  communications,  space  is  said  to  be  the  final  frontier.  It  has 
been  shown  that  MIMO  systems  and  techniques  can  lead  to  a  big  increase  in 
system  performance  in  terms  of  capacity,  coverage  and  signal  quality,  all  of 
which  ultimately  lead  to  increased  spectrum  efficiency.  While  the  future  is 
promising,  there  are  still  problems  that  need  to  be  solved.  The  goal  of  this 
report  is  to  propose  a  new  class  of  turbo  space- time  codes  with  the  application 
of  TVLT.  Compared  to  the  existing  space-time  codes,  the  new  codes  proposed 
in  this  report  can  be  designed  with  a  significantly  lower  complexity  procedure, 
while  achieving  a  higher  data  rate  and  better  performance.  The  new  codes 
can  also  be  extended  and  applied  to  systems  with  a  very  high  data  rate  and 
a  large  number  of  transmit  antennas. 

The  major  contributions  of  this  report  include  the  following. 

•  The  effects  of  the  TVLT  technique  on  the  diversity  gain  and  coding 
gain  of  turbo  space-time  codes  are  studied.  It  has  been  shown  that 
with  TVLT  full  diversity  and  good  coding  gain  is  achieved  with  high 
probability.  The  probability  that  these  design  goals  are  not  met  de¬ 
creases  exponentially  with  the  minimum  Hamming  distance  between 
the  codewords.  Hence,  the  need  to  test  every  pair  of  codewords  for  the 
rank  criterion  and  coding  gain  criterion  [1]  is  eliminated.  The  design 
criterion  for  turbo-STC  with  TVLT  is  to  maximize  the  Hamming  dis¬ 
tance.  Turbo-STCs  are  thus  good  candidates  to  be  used  with  TVLT 
to  achieve  good  coding  gain,  since  they  usually  have  large  minimum 
Hamming  distance  following  the  existing  design  guidelines.  Hence,  the 
design  of  full  diversity,  good  coding  gain  turbo-STCs  is  significantly 
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simplified. 


•  The  effect  of  TVLT  on  the  achievable  rate  of  turbo-STC  is  studied. 
It  is  shown  that  while  achieving  full  diversity,  turbo-STC  can  achieve 
much  higher  data  rate  than  what  can  be  achieved  by  any  full  diversity 
space-time  codes  with  constant  modulation  [1],  A  new  rate  limit  for 
turbo-STC  with  TVLT  is  proposed. 

•  The  convenient  and  powerful  EXIT  chart  method  [13]  is  extended  to 
analyze  and  predict  the  performance  of  turbo-STC  with  TVLT.  It  is 
shown  that  with  TVLT  the  performance  of  turbo-STC  is  only  deter¬ 
mined  by  the  singular  values  of  the  channel  gain  matrix.  This  reduction 
in  effective  channel  parameters  enables  a  low  dimension  description  of 
the  MIMO  channel  conditions,  which  further  enables  a  low  complexity 
search,  with  the  EXIT  chart  method,  for  the  channel  conditions  over 
which  the  turbo-STC  decoder  will  converge.  The  FER  performance  of 
turbo-STC,  can  then  be  predicted  with  significantly  lower  complexity 
than  which  is  needed  for  a  full  simulation.  The  low  complexity  per¬ 
formance  analysis  and  prediction  can  significantly  expedite  the  code 
comparison  and  selection  in  turbo-STC  design. 

•  Multilevel  turbo-STC  with  TVLT  is  proposed  to  increase  the  data  rate 
without  exponentially  increasing  the  decoding  complexity.  Multidimen¬ 
sional  set  partition  is  applied  which  enables  higher  level  codes  to  have 
higher  data  rate  and/or  lower  trellis  complexity.  An  iterative  multi¬ 
stage  decoder  is  proposed  and  reduced  complexity  decoder  with  a  list 
sphere  decoder  is  also  proposed  as  an  alternative. 

•  Layered  turbo  space-time  (LTST)  coding  is  proposed  for  MIMO  sys¬ 
tems  with  large  number  of  transmit  antennas,  to  avoid  exponentially 
increasing  decoding  complexity  with  the  number  of  transmit  antennas. 
Two  different  structures,  V-LTST  and  H-LTST,  are  proposed  for  sys¬ 
tems  with  low  and  high  order  modulation,  respectively.  In  V-LTST  a 
very  low  complexity  interference  cancellation  technique  is  used  while 
superior  performance  is  achieved.  In  H-LTST,  ordered  SIC  is  imple¬ 
mented,  enabled  by  the  TVLT  application  to  each  layer.  Significant 
performance  improvement  is  demonstrated  compared  to  the  conven¬ 
tional  layered  space-time  coding  structure  where  ordering  is  not  avail¬ 
able  in  the  SIC  detection. 
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A  group  of  new  codes  is  constructed  based  on  the  design  techniques 
proposed  in  this  report.  Their  simulation  results  are  reported. 


Although  much  study  has  been  done  in  this  report  on  turbo-STC  with 
TVLT  scheme,  the  potential  of  this  scheme  are  far  from  being  fully  explored. 
Anticipated  future  work  may  contain  the  following. 

•  Use  turbo-STC  with  TVLT  scheme  to  construct  adaptive  modulation 
and  coding  (AMC)  MIMO  systems.  With  TVLT,  whether  the  turbo- 
STC  decoder  will  converge  or  not  can  be  accurately  predicted  based  on 
the  channel  state  information.  The  receiver  can  inform  the  transmitter 
to  select  a  code  with  the  maximum  rate  that  guarantees  the  convergence 
of  the  decoder. 

•  Evaluate  the  possibility  of  applying  turbo-STC  with  TVLT  as  a  mul¬ 
tiple  access  method  for  multi-user,  MIMO  systems.  The  sequence  of 
pseudo-random  unitary  matrices  used  in  TVLT  can  be  regarded  as  the 
random  codes  for  different  users  in  CDMA  systems. 

•  Study  the  applications  of  turbo-STC  with  TVLT  to  construct  secured 
MIMO  wireless  systems.  With  the  TVLT,  the  transmitted  signals  are 
multiplied  by  time  varying  pseudo-random  unitary  matrices  and  they 
appear  random.  The  signals  will  be  difficult  to  decode  unless  the  ran¬ 
dom  matrices  are  known. 

•  Evaluate  the  performance  of  turbo-STC  with  TVLT  in  multiple  input 
multiple  output  orthogonal  frequency  division  multiplexing  (MIMO- 
OFDM)  systems.  OFDM  technology  is  playing  a  more  and  more  im¬ 
portant  role  in  broadband  wireless  communications.  It  has  been  ap¬ 
plied  to  standards  IEEE  802.11a/g/n  and  DVB-H,  and  will  be  used  in 
other  standards  such  as  IEEE  802.16  and  802.20.  Over  the  frequency 
selective  channel,  the  minimum  Hamming  distance  requirement  of  the 
code  used  with  TVLT  may  be  lowered.  If  this  is  true,  the  application 
of  TVLT  can  be  extended  to  conventional  STC  coded  MIMO-OFDM 
systems. 
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